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The  model  presented  in  this  dissertation  incorporates  strong  collision  effects  in  the 
electron  broadening  of  spectral  lines  from  charged  radiators  in  hot,  dense  plasmas.  This 
model  is  an  extension  of  an  existing  model,  often  referred  to  as  the  unified  theory,  that 
was  originally  developed  to  calculate  frequency-dependent  electron-broadening  effects 
consistently  from  spectral  line  center  to  the  far  line  wings  for  neutral  hydrogen  radiators 
and  later  adjusted  for  singly-charged  helium  radiators.  Specifically,  the  extension  allows 
for  the  consistent  calculation  of  electron-broadening  effects  on  line  spectra  emitted  from 
mid-Z  radiators  (i.e. , argon).  This  is  accomplished  through  the  use  of  relativistic  radiator 
state  data  and  also  the  inclusion  of  strong  collision  effects  by  incorporating,  in  addition 
to  the  dipole  contribution,  the  effects  of  the  monopole  term  and  all  other  contributing 
higher-order  “multipole”  terms  in  the  electron-radiator  interaction.  This  model  preserves 
the  semiclassical  treatment  of  the  original  model  where  the  perturbing  electrons  move  in 
classical  trajectories  and  the  radiator  is  treated  through  quantum  mechanical  means.  It 
also  preserves  the  calculation  of  the  electron-broadening  operator  to  all-orders  in  terms 
of  the  electron-radiator  interaction. 
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Line  spectra  using  the  full-Coulomb,  semiclassical,  all-order  electron-broadening 
model  are  compared  to  line  spectra  from  a dipole-approximation  version  of  this 
model.  These  line  spectra  are  also  compared  to  line  spectra  from  full-Coulomb  and 
dipole-approximation  versions  of  a full  quantum  mechanical  electron-broadening  model 
expanded  to  second-order  in  a Born-like  expansion  in  terms  of  the  electron-radiator 
interaction.  Significant  shifts  and  distortions  in  the  full-Coulomb  line-shape  calculations 
as  compared  to  the  dipole-approximation  line-shape  calculations  at  electron  densities 
greater  than  1 x 1024cm-3  show  that  a dipole  approximation  of  the  electron-radiator 
interaction  is  no  longer  reasonable.  Differences  in  the  amount  of  shift  and  distortion 
in  the  comparisons  of  full-Coulomb  calculations  between  the  semiclassical,  all-order  and 
quantum  mechanical,  second-order  models,  specifically  in  the  calculation  of  merging 
spectral  lines  at  electron  densities  above  1 x 1024  cm~3,  suggest  that  an  all-order  model 
is  needed  as  the  density  increases. 
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CHAPTER  1 
INTRODUCTION 


When  a radiating  ion  is  immersed  in  a dense  plasma,  the  energy  levels  of  its  bound 
states  are  modified  by  interactions  with  the  surrounding  ions  and  electrons.  This  results 
in  spectral  lines  that  are  significantly  broadened  and  shifted  from  the  spectral  lines 
of  an  isolated  radiator  [1,2],  The  magnitude  of  this  effect  depends  on  the  identity  of 
the  radiator  and  the  density  and  temperature  of  the  ions  and  electrons.  Therefore, 
the  width  and  shift  of  the  lines  provide  a useful  diagnostic  measurement  of  the  density 
and  temperature  of  hot,  dense  plasmas  [3-5].  As  the  density  of  a plasma  increases, 
models  describing  the  broadening  and  shifting  of  the  spectra  must  be  modified  to  more 
accurately  account  for  high-density  effects.  This  dissertation  specifically  examines  the 
effects  of  strong  electron  collisions  on  the  spectra  of  a radiator  immersed  in  a plasma 
of  ions  and  electrons  and  the  potential  use  of  the  resulting  spectra  as  a density  and 
temperature  diagnostic  in  the  analysis  of  experimental  data.  A special  focus  is  placed 
on  plasma  conditions  in  which  some  of  the  spectral  lines  begin  to  merge  together. 

It  is  appropriate  at  this  point  to  parse  the  title  of  this  work,  as  it  contains  some  terms 
that  need  to  be  explained.  In  the  context  of  this  work,  a “strong  collision”  refers  to  those 
collisions  in  which  the  strength  of  the  interaction  between  the  radiator  and  a perturbing 
electron  is  large  enough  and  occurs  over  a sufficient  period  of  time  to  disrupt  a radiative 
transition  from  an  upper  state  to  a lower  state.  This  work  does  not  include  the  effects  of 
ionizing  collisions  and  electron  capture.  “Electron  broadening”  refers  to  the  broadening 
effects  of  perturbing  electrons  through  their  modifications  to  the  bound  states  of  the 
radiator.  “Charged  radiators”  corresponds  to  ions  in  the  range  of  6 < Z < 36  that  are 
undergoing  radiative  transitions  and  emitting  spectra.  “Hot,  dense  plasmas”  refers  to 
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nondegenerate  plasmas  where  density  broadening  effects,  specifically  Stark  broadening 
and  electron  broadening,  are  the  dominant  broadening  mechanisms  on  the  spectra  and 
are  large  enough  to  cause  spectral  lines  to  overlap,  or  “merge”,  together. 

Improvements  in  laser  technology  over  the  past  20  years  led  to  the  creation  of  plas- 
mas with  increasingly  higher  densities  and  temperatures  [6-8].  To  analyze  and  under- 
stand these  plasmas,  original  theoretical  line  broadening  models  have  been  advanced 
and  improved  upon.  Early  models  began  with  the  Stark  broadening  of  spectra  from 
neutral  radiators  in  the  presence  of  an  electric  field  from  the  distribution  of  the  sur- 
rounding ions,  commonly  referred  to  as  the  ion  microfield  [9,10],  along  with  an  electron- 
broadening model  utilizing  a dipole  interaction  approximation.  As  the  density  and 
temperature  of  the  plasmas  increased,  the  study  of  spectra  from  ionic  radiators  became 
necessary  [11-13]  because  the  hydrogen  atoms  became  completely  ionized,  thus  requiring 
a multipole,  or  full-Coulomb  calculation  of  the  interaction  between  the  charged  radiator 
and  the  perturbing  plasma  electrons  in  the  electron-broadening  calculation  [14,15].  The 
inclusion  of  the  monopole  term  in  the  interaction  led  to  the  prediction  of  an  experi- 
mentally noticeable  red  shift  in  the  location  of  the  spectral  lines  as  the  density  of  the 
plasmas  increased  [16-19].  This  shift  was  shown  to  increase  with  increasing  principal 
quantum  number  and  electron  density  and  to  decrease  with  increasing  temperature. 

As  a contextual  aid  for  the  model  presented  in  this  dissertation,  we  begin  by  briefly 
describing  the  experimental  program  currently  underway  to  generate  the  experimental 
line  spectra  that  are  used  to  probe  the  plasma  conditions  achieved  in  the  experiment 
through  comparisons  with  theoretical  line  spectra  from  different  models.  In  order  to 
generate  plasmas  with  the  high  electron  densities  needed  to  observe  high-density  effects, 
a high-powered  laser,  specifically  the  64  beam  OMEGA  laser  system  at  the  Laboratory 
for  Laser  Energetics  at  the  University  of  Rochester,  is  used  to  implode  a small  plastic 
shell  (CH  microballoon)  filled  with  a gas  consisting  of  deuterium  with  trace  amounts  of 
a mid-Z  element  (i.e.,  argon).  The  main  variables  of  interest  in  our  experiments  include 
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the  shell  thickness  of  the  CH  microballoon,  the  microballoon’s  diameter,  the  pressure 
and  composition  of  the  fill  gas,  the  laser  beam  pulse  shape,  and  the  total  energy  delivered 
on  target. 

For  our  most  recently  analyzed  experiments  through  our  NLUF  (National  Laser  User 
Facility)  campaign,  the  CH  shell,  or  ablator,  thickness  generally  was  20  microns,  and 
the  CH  inner  shell  diameter  was  940  microns.  The  composition  and  pressure  of  the 
fill  gas  ranged  from  0.3%  to  2%  atomic  fraction  argon  in  7-20  Atm  of  D2.  The  laser- 
beam  pulse  shape  was  either  a 1 ns  ramp  to  1 ns  flat-top  pulse  or  a 1 ns  flat-top  pulse. 
The  total  ultraviolet  laser  energy  delivered  on  target  was  in  the  range  of  19-24  kJ.  The 
energy  delivered  on  target  would  generally  be  a bit  larger,  but  some  total  beam  energy 
was  sacrificed  to  achieve  excellent  beam  balance  of  about  A E/E  ~ 3%  RMS.  The 
experiments,  with  a few  exceptions,  were  successful.  The  framing  camera  images  and 
the  static  pinhole  images  all  showed  nearly  circular  core  images,  showing  no  significant 
excursions  from  a spherical  implosion.  Good  argon  K-shell  data  were  observed  in  most 
of  the  experimental  shots. 

Our  primary  diagnostic  was  a streaked  x-ray  crystal  spectrometer,  or  SSCA,  operat- 
ing with  a flat  RbAP  (rubidium  acid  phthalate)  crystal  and  an  Au  photocathode.  The 
instrument  provided  approximately  a 1.9  ns  window  with  a resolution  of  25  ps.  The 
spectra  were  recorded  on  film  and  then  digitized.  Because  the  isolated  line  position 
is  very  important  in  determining  important  dense-plasma  effects,  including  line  shifts 
and  line-merging  effects,  corrections  for  film  sensitivity,  filtering,  photocathode  response, 
streak-camera  distortions  including  curvature  of  isotemporal  lines  and  streak  angle,  and 
crystal/photocathode  irregularities  were  performed  on  the  data.  Lineouts  averaging  over 
25  ps  were  then  extracted  from  the  images,  and  the  resulting  time-dependent  spectra 
exhibited  time-varying  Ar  K-shell  data.  Spectral  dispersion  was  determined  using  early- 
time  lineouts  when  the  electron  density  was  small. 
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Our  analysis  of  the  data  proceeds  by  performing  a least  squares  (y2)  fit  to  the  spectral 
data  using  our  group’s  fitting  model.  The  model  contains  Stark  broadened  line  shapes, 
NLTE  population  distributions,  and  corrections  for  the  transfer  of  the  thick  alpha  lines 
and  the  opacity  of  the  diagnostic  lines.  The  free  parameters  in  our  fitting  model  are  the 
emissivity-averaged  core  electron  temperatures  and  densities.  This  fitting  process  leads 
to  an  inference  of  time-dependent  core  conditions,  something  useful  in  its  own  right. 

The  fitting  analysis  of  our  most  recently  analyzed  experiments  confirmed  that  the 
experimental  shots  provided  plasma  conditions  approaching  those  needed  for  our  investi- 
gation of  dense-plasma  effects.  We  regularly  reached  densities  just  above  2 x 1024  cm-3. 
At  this  density,  high-density  deviations  from  a linear  Stark-broadening  theory  that  does 
not  include  mixing  between  the  upper-state  manifolds  of  adjacent  spectral  series  mem- 
bers, specifically  in  line  merging  behavior  and  strong-collision  effects,  are  subtle.  How- 
ever, we  saw  indications  in  the  data  that  a nonlinear  Stark-broadening  theory  that 
includes  mixing  between  the  upper-state  manifolds  of  adjacent  series  members  yielded 
fits  with  a somewhat  lower  y2  than  those  from  the  linear  theory,  especially  in  the  region 
of  the  Ly-y  line,  as  is  shown  in  Figure  1.1.  As  our  NLUF  campaign  progresses,  we  are 
confident  that  changes  in  the  shell  thickness,  fill  gas  pressure,  and  pulse  shape  should 
give  densities  on  the  order  of  4 x 1024  cm-3  where  the  differences  are  not  so  subtle. 

As  was  just  noted,  we  can  currently  generate  electron  densities  of  a magnitude  such 
that  strong  collisions  are  beginning  to  affect  the  spectra  from  hydrogenic  and  helium- 
like argon.  In  other  words,  the  spectra  become  sufficiently  broad  and  shift  to  the  extent 
that  the  contributions  from  strong  collisions  must  be  studied  in  greater  detail.  Previous 
models  that  limit  the  effects  of  strong  collisions  by  including  only  those  contributions 
to  second-order  in  the  perturbing  electron-radiator  interaction  are  valid  over  a range 
of  A u)  < Zujpiasma,  where  A uj  is  the  frequency  separation  from  line  center  [20].  For 
example,  with  current  electron  densities,  the  Lyman-5  line  of  Argon  broadens  and  shifts 
to  the  extent  that  it  has  already  merged  with  the  Lyman-q  line.  In  fact,  the  Lyman-q5 
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strates  that  a quantum-mechanical  nonlinear  model  to  second-order  in  the  full-Coulomb 
electron-radiator  interaction  gives  a slightly  better  fit  than  the  linear  version  of  this 
model  at  an  electron  density  of  1.9  x 1024  cm-3  and  a temperature  of  1.15  keV.  The 
value  of  x2  is  slightly  smaller  for  the  nonlinear  fit. 

spectral  feature  is  already  affecting  the  Lyman-/!  line  significantly.  In  addition,  these 
merged  spectral  features  tend  to  be  significantly  broader  than  2upiasma . 

To  incorporate  strong  collisions,  as  defined  above,  into  the  theoretical  foundation 
of  plasma  line  broadening  for  mid-Z  ionic  radiators,  we  develop  an  electron-broadening 
model  similar  to  a model  developed  by  Vidal,  Cooper,  and  Smith  (VCS)  back  in  the 
late  1960s  and  early  1970s  for  neutral  radiators  [21-24],  and  extended  a few  years  later 
by  Greene  and  Cooper  for  ionized  helium  [25-27].  The  VCS  model  was  originally  devel- 
oped to  calculate  spectral  lines  of  neutral  hydrogen  radiators  from  line  center  all  the 
way  out  to  the  far  line  wing  consistently.  This  model  is  sometimes  referred  to  as  the 
unified  theory  of  line  broadening  because  it  consistently  joined  the  impact  model  of 
line  broadening  that  is  appropriate  around  line  center  to  the  one-electron  model  of  line 
broadening  that  is  appropriate  in  the  far  wings  of  a line.  The  VCS  model  is  a frequency- 
dependent,  binary-collision  electron-broadening  model  that  includes  terms  to  all-order 
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in  the  interaction  between  the  radiator  and  the  perturbing  electron  through  the  calcula- 
tion of  a time  development  operator.  It  is  a semiclassical  model,  meaning  that  it  treats 
the  perturbing  electrons  as  classical  particles  using  classical  trajectories  while  treating 
the  radiator  through  quantum  mechanical  means.  Because  the  VCS  model  focuses  on 
low-Z  neutral  radiators  and  excludes  strong-collision  contributions,  this  model  uses  only 
the  dipole  term  in  the  perturbing  electron-radiator  interaction  and  uses  nonrelativistic 
wavefunctions  to  describe  the  radiator  states.  A no-quenching  approximation  is  also 
applied  in  the  VCS  model  so  that  perturbations  on  the  radiator  do  not  induce  non- 
radiative  transitions  between  the  initial  and  final  states.  Electron  correlations  in  this 
model  are  included  by  considering  only  those  perturbing  electron-radiator  interactions 
which  occur  within  a sphere  about  the  radiator  with  a radius  of  the  Debye  length.  Any 
electrons  outside  of  this  sphere  do  not  interact  with  the  radiator  because  the  radiator  is 
effectively  screened  by  the  electrons  within  the  Debye  sphere. 

The  difference  between  the  VCS  model  and  the  model  as  modified  by  Greene  and 
Cooper  was  to  replace  the  straight-line  trajectories  of  the  perturbing  electrons  for  neutral 
radiators  with  hyperbolic  trajectories  for  charged  radiators.  The  other  aspects  of  the 
model  basically  remained  unchanged.  It  is  this  version  of  the  VCS  model  modified  for 
charged  low-Z  radiators  that  is  similar  to  the  model  developed  in  this  dissertation. 

In  order  to  calculate  theoretical  spectra  for  higher-Z  radiators  (i.e.,  Argon)  and  to 
include  contributions  from  strong  collisions,  the  model  presented  in  this  work  differs 
from  the  VCS  charged-radiator  model  in  these  important  ways: 

• The  interaction  between  the  perturbing  electron  and  the  radiator  is  calculated 
using  a multipole,  or  full-Coulomb  calculation. 

• Relativistic  wavefunctions  and  the  corresponding  energy  level  structure  from  a 
relativistic  atomic  physics  package  are  employed  to  describe  the  atomic  physics  of 
mid-Z  radiators. 

• A minimum  separation  cutoff  is  incorporated  into  the  model.  This  cutoff  is  the 
Landau  length  modified  by  the  radiator  charge  and  quantum  diffraction  effects. 
The  cutoff  corrects  for  the  singularity  caused  by  the  treatment  of  the  electrons  as 
classical  particles. 
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• The  distribution  function  is  generalized  to  provide  temperature-dependent  changes 
in  the  relative  populations  of  the  radiator  states  modified  by  the  interaction 
between  the  perturbing  electron  and  the  radiator. 

• The  no-quenching  approximation,  which  prohibits  nonradiative  transitions  caused 
by  perturbations  between  states  of  differing  principle  quantum  number  n,  must  be 
relaxed  to  allow  for  nonradiative  transitions  between  upper  states  of  spectral  lines 
which  are  beginning  to  merge  together. 

With  these  improvements,  this  model  is  able  to  accurately  calculate  isolated  and  merging 
spectral  lines  from  mid-Z  radiators. 

In  addition  to  the  calculation  of  isolated  and  merging  spectral  lines,  there  are  several 
other  applications  for  this  model.  Differences  in  merged  spectral  lines  at  higher  electron 
densities  due  to  strong-collision  effects  from  the  higher-order  terms  has  the  possibility 
of  altering  what  we  perceive  as  the  distinction  between  continuum  and  bound  states. 
Because  ionization  balance  is  sensitive  to  any  changes  in  the  energy  level  structure  [28-30] 
and  opacity  is  sensitive  to  the  energy  location  and  intensity  of  line  emission,  an  under- 
standing of  the  effects  of  strong  collisions  on  the  merged  composite  spectral  features  that 
lie  near  the  continuum  edge  can  be  important. 

Even  at  lower  densities,  the  wings  of  a line  profile  are  significantly  dependent  on 
the  effects  of  strong  collisions  through  the  higher-order  terms.  This  can  have  significant 
effects  on  Rosseland  group  opacities  [31].  They  are  very  sensitive  to  changes  in  the 
intensity  of  the  line  profile  in  the  wings  because  they  are  proportional  to  a harmonic 
mean  of  the  absorption  coefficient.  Rosseland  group  opacities,  in  turn,  are  important  in 
the  calculation  of  multi-group  diffusion  radiative  transfer  effects.  These  effects,  which  are 
important  in  determining  the  amount  of  radiation  that  passes  through  different  portions 
of  the  frqequency  spectrum,  can  be  of  importance  in  both  astrophysical  studies  and  in 
fusion  research. 

To  outline  this  dissertation,  in  Chapter  2,  we  present  a model  of  spectral  line  broad- 
ening for  radiators  immersed  in  hot,  dense  plasmas.  We  consider  several  broadening 
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• Doppler  broadening  due  to  thermal  motion  of  the  radiator  is  included. 

• Stark  broadening  due  to  the  distribution  of  static-ion  fields  is  included. 

• Dynamical  ion  broadening  due  to  transitions  between  the  static-ion  fields  is 
included. 

• Electron  broadening  due  to  strong  electron-radiator  collisions  is  included. 

For  the  plasma  conditions  under  consideration,  the  dominate  broadening  mechanisms 
are  the  Stark  broadening  due  to  the  distribution  of  static-ion  fields  and  the  electron 
broadening  due  to  electron-radiator  collisions.  The  relative  importance  of  these  broad- 
ening mechanisms  are  discussed  in  detail  by  others  [18],  and  is  not  revisited  here.  Some 
simplifying  assumptions  are  made  for  calculational  reasons  which  will  need  to  be  further 
investigated  in  the  future. 

Chapter  3 focuses  on  the  calculation  of  the  electron-broadening,  or  width-and-shift, 
operator.  We  introduce  a convenient  form  for  the  electron-radiator  interaction,  discuss 
the  classical-path  approximation,  and  apply  some  approximations  to  allow  for  tractable 
calculations  of  this  operator. 

In  Chapter  4,  we  discuss  the  details  of  the  calculation  of  spectral  features.  We  then 
provide  an  analysis  of  widths  and  shifts  due  to  electron  broadening.  We  also  show 
comparisons  of  line  shapes  using  the  electron-broadening  model  described  in  Chapter 
3 with  line  shapes  utilizing  different  electron-broadening  models.  Finally,  we  look  at 
comparisons  of  results  from  calculations  involving  the  merging  of  spectral  lines. 

In  the  Conclusions,  we  summarize  our  findings  and  consider  the  effects  of  the  treat- 
ment of  strong  electron  collisions  in  our  model  on  other  areas  in  the  field  of  dense  plasma 
physics. 


CHAPTER  2 

LINE  BROADENING  IN  HOT  DENSE  PLASMAS 


The  first  step  in  considering  line  broadening  in  hot,  dense  plasmas  is  to  develop 
a model  describing  the  plasma.  This  model  is  crucial  to  an  understanding  of  how  the 
interactions  of  the  plasma  with  a highly  ionized  radiator  affects  the  spectral  lines  emitted 
by  the  radiator. 

The  radiators  (i.e.,  argon)  are  assumed  to  be  well  separated  in  the  plasma  and  do 
not  interact  with  each  other.  Thus,  the  plasma  is  broken  up  into  noninteracting  cells 
where  each  cell  contains  one  radiator  along  with  electrons  and  ions  that  perturb  the 
quantum  states  of  the  radiator.  It  is  from  within  this  perturbing  plasma  environment 
that  the  radiator  emits  its  radiation.  The  frequency-dependent  power  distribution  of  the 
radiation,  specifically  electric  dipole  radiation,  can  be  written  in  terms  of  an  ensemble 
average  over  all  the  possible  states  of  radiator  together  with  the  plasma  environment 
as  [1,2] 

A.  .4 

P M = i^'52\{b\de-'k-r\a)\2  pa5(u  - uab)  (2.1) 

a,  b 

where  a represents  an  initial  state  of  the  radiator  and  plasma  environment,  weighted  by 
the  corresponding  eigenvalue  of  the  density  matrix  pa , b represents  a final  state  of  the 
radiator  and  plasma  environment,  and  hujab  is  the  energy  of  the  transition.  The  electric 
dipole  operator  of  the  system  is  denoted  by  d , the  vector  k represents  the  wave  vector 
of  the  emitted  radiation,  and  r is  the  center  of  mass  position  of  the  radiator. 

Because  to4  does  not  vary  significantly  over  the  frequency  range  of  the  spectra  of 
interest,  the  line  shape  is  defined  from  the  power  distribution  formula  as 

I(u)  = ]£<a|  eikrd\  b)  ■ {b\de~ikr\ a)pa5{uj  - uab) . (2.2) 

a,  b 
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Using  the  definitions 

hu>ab  = Ea-Eb  and  dk  = de~'kr  (2.3) 

and  the  integral  representation  of  the  Dirac  delta  function,  the  line  shape  can  be  rewrit- 
ten in  the  form 

I(u)  = — / dt  ei(“~i£)t  Tr{dfc  • e-im/hpd]k  eiHt/n}  (2.4) 

27r  J- oo 

where  the  trace  is  over  the  plasma-radiator  system.  The  expression  inside  the  trace 
is  the  autocorrelation  function  of  the  radiator  dipole,  keeping  in  mind  that  the  trace  is 
invariant  over  cyclic  permutations.  The  damping  factor  e in  the  transform  represents  the 
natural  line  width  which  corresponds  to  the  lifetime  of  the  upper  state  of  the  transition. 
The  density  matrix  operator  has  the  form 

_ exp  (~H/kBT) 

Tr{exp (-H/kBT)Y 

To  complete  the  definition  of  the  line  shape  function,  the  form  of  the  Hamiltonian 
H must  be  specified.  Separating  out  the  important  interactions  in  the  plasma- radiator 
system,  H has  the  form 

H = Hrcm  + Hrint  + Hi  + He  + Vi>e  + VI, r + Ue,r  . (2.6) 

The  first  term  refers  to  the  center-of-mass  motion  of  the  radiator  while  the  second  term 
refers  to  the  internal  states  of  the  radiator.  The  next  two  terms  refer  to  the  sum  of  the 
kinetic  energies  and  the  interactions  within  the  respective  subsystem,  or 

Hi  = Ki  + Vij  and  He  = Ke  + Ve>e  . 


(2.7) 
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The  last  three  terms  represent  the  interactions  between  the  radiator,  the  electrons,  and 
the  ions. 

Now  that  we  have  the  starting  expression  for  the  line  shape,  we  can  make  certain 
approximations  based  on  the  properties  of  the  plasmas  considered.  Because  these  hot 
plasmas  have  a radiator  (i.e. , argon)  that  is  usually  much  more  massive  then  the  per- 
turbing ions  (i.e.,  deuterium)  and  obviously  much  more  massive  then  the  electrons,  it  is 
usually  appropriate  to  assume  that  the  center-of-mass  motion  of  the  radiator  is  statis- 
tically independent  from  the  other  degrees  of  freedom  [32-34] . Stated  another  way,  the 
impulse  imparted  to  the  radiator  from  the  electron  and  ion  perturbations  is  small  rela- 
tive to  the  magnitude  of  the  radiator’s  momentum  over  the  lifetime  of  the  upper  state 
of  the  transition.  Assuming  that  the  radiators  are  well  separated  and  can  be  treated 
as  an  ideal  gas,  the  line  shape  can  be  expressed  as  a convolution  of  a Doppler  profile 
(gaussian)  and  a Stark- broadened  line  shape  in  the  radiator  center-of-mass  frame  [35,36] 


where 


2cr ' 


a/2  (2kBT 


mr 


(2.9) 


The  line  shape  in  the  radiator  center-of-mass  frame  is  given  by 


i r °°  , , 

ICM(u)  = — / d t ei(u~ie)t  Tr{d  • e~iH  t^np'dieiH  t'h]  (2.10) 

2tt  J _«*, 


where  the  density  matrix  operator  now  has  the  form 

exp(—H/kBT) 
Tr{exp  (-H'/kBT)} 


P 


(2.11) 
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with  the  Hamiltonian  given  by 


(2.12) 


To  keep  the  notation  simple,  we  drop  the  internal  state  indicator  on  the  isolated  radiator 
Hamiltonian. 


From  this  point  on,  the  center-of-mass  notation  on  the  Stark-broadened  line  shape  func- 
tion is  suppressed.  The  notation  I (to)  represents  the  Stark-broadened  center-of-mass 
line  shape. 

As  a lead  in  to  the  simplification  of  the  line  shape  expression,  it  is  convenient  to 
regroup  the  terms  of  the  Hamiltonian  in  the  form 


The  center-of-mass  line  shape  can  be  simplified  by  taking  advantage  of  the  fact  that 
the  expression  inside  the  trace  is  the  autocorrelation  function  of  the  electric  radiator 
dipole,  C(t).  Because  d is  hermitian,  C(t)  has  the  property 


c(-f)  = C(ty 


(2.13) 


that  can  be  employed  to  write  the  line  shape  as 


1 r °°  , , 

ICmH  = -ft  / d t e'^-'t)trTxrJd  ■ e~iH  t'hp'deiH  t/h)  . (2.14) 

7T  Jn 


Jo 


(2-15) 


where  the  ion  Hamiltonian  is  given  by 


(2.16) 
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and  the  electron  Hamiltonian  is  given  by 


(2.17) 


Note  that  the  monopole,  or  point-charge,  interactions  between  the  radiator  and  the 
perturbing  ions  and  electrons,  which  are  denoted  by  the  superscript  “ 0 ”,  have  been 
separated  out  and  combined  with  the  ion  and  electron  portions  of  the  Hamiltonian, 
respectively.  Thus,  the  dynamical  behavior  of  the  ions  and  electrons  in  the  given  radiator 
cell  of  the  plasma  is  defined  in  the  presence  of  the  charge  of  the  radiator.  Due  to  the 
repulsion  between  the  charged  radiator  and  the  perturbing  ions,  it  should  be  a reasonable 
approximation  that  the  ions  do  not  penetrate  the  the  orbitals  of  the  radiator.  The 
remainder  of  the  ion-radiator  interaction  can  be  approximated  by  a radiator  dipole 
interaction  with  the  electric  field  due  to  the  ions 


interaction.  In  other  words,  this  term  represents  the  interaction  of  the  perturbing  elec- 
trons with  the  internal  electronic  structure  of  the  radiator. 

Even  without  correlations  between  Doppler  broadening  and  Stark  broadening  in  the 
line  shape  expression,  the  calculation  of  the  Stark-broadened  center-of-mass  line  shape 
function  is  still  very  complicated.  A method  to  simplify  the  calculation  of  the  line  shape 
function  is  employed  by  integrating  the  line  shape  at  a specific  value  of  the  static-ion 


(2.18) 


The  effects  of  the  higher-order  terms  of  on  various  line  shapes  have  been  evalu- 
ated [37]  and  appear  to  have  a relatively  small  effect  in  comparison  to  the  effects  of 
plasma-induced  line  shifts  on  the  line  shapes.  The  term  V$  now  consists  of  the  interac- 
tion of  the  radiator  with  the  perturbing  electrons  excluding  the  point-charge  monopole 
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field  over  a distribution  of  static-ion  fields 

I(uj)  = J d eQ(e)J(u>,e)  (2-19) 

where 

Q(e)  = Trt{Pi5(e  - Esi )}  (2.20) 

is  the  probability  that  the  distribution  of  ions  around  the  radiator  results  in  an  electric 
field  (microfield)  of  the  value  e at  the  location  of  the  radiator.  The  functional  form  of 
the  probability  distribution  function  is  shown  in  Figure  2.1  for  various  temperatures  and 
densities.  The  electron-broadened  line  shape  evaluated  at  a given  value  of  the  static-ion 
microfield,  is  given  by 


J(u,e)  = --3Trr{d-  D{co,e)}  (2.21) 

7T 

and  its  functional  form  is  shown  in  Figure  2.2.  The  details  of  the  computation  of  Q(e)  and 
are  described  at  the  beginning  of  Chapter  4.  A comparison  of  these  expressions 
to  equation  2.14  shows  that  D(lj,  e)  has  the  form 

D(w,  e)  = -f-  / d t Trp  { 5 (e  - Esi)  p' e~'H  t/nd  eiH  t/h  \ (2.22) 

Q\e)  Jo  L 

where  the  distribution  function  of  the  plasma-radiator  system  p commutes  with  the 
time  evolution  operator  of  the  plasma-radiator  system  e~lH 

The  next  step  in  the  calculation  of  the  line  shape  expression  involves  the  derivation 
of  the  equation  of  motion  for  D(u,  e)  along  with  the  derivation  of  its  solution.  The 
details  of  this  derivation  are  given  in  Appendix  A,  and  the  resulting  equation  of  motion 
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Figure  2.1:  Microfield  probability  distribution  functions  47re2Q(e)  are  shown  for  1% 
atomic  fraction  argon  in  D2.  The  microfield  is  given  in  units  of  e0  = e2/r0  where  r0  is  the 
inter-electron  spacing,  a)  The  density  behavior  of  47re2Q(e)  is  shown  at  a temperature 
of  1 keV.  Note  that  the  most  probable  field  value  increases  with  ne  keeping  in  mind  the 
density  dependence  of  the  reduced  units,  b)  The  temperature  behavior  of  47re2<5(e)  is 
shown  at  a density  of  1 x 1024  cm-3.  Note  that  the  most  probable  field  value  increase 
with  temperature  as  well. 
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Figure  2.2:  J(oj,  e)  is  shown  for  the  Lyman-a:  line  of  Ar+17  at  a density  of  1 x 1024  cm-3 
and  a temperature  of  1 keV.  Field  values  are  given  in  units  of  e0  = e2/ro  where  ro  is  the 
inter-electron  spacing,  a)  This  profile  is  calculated  using  the  field  value  e = 0.  Thus, 
there  is  no  ion  Stark  splitting,  b)  This  profile  is  calculated  using  the  field  value  e = 1.1, 
the  most  probable  field  for  these  plasma  conditions. 
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has  the  form 


u-Lr-  L>V(c)  - Be(e)  - Me(cu,  e)j  D(u,  e) 

- J de' Mi(e,e'-,uj)D(u>,e')  = p€iTd  (2.23) 


where  we  have  assumed  that  there  are  no  dynamical  correlations  between  the  electrons 
and  the  ions.  The  broadening  contribution  from  the  isolated  radiator  operator  Lr  is 
given  by 

LrD(u,e)  = ^[Hr,D(u,e)}  (2.24) 

and  the  broadening  contribution  from  the  partial  removal  of  the  degeneracy  of  the 
radiator  states  in  the  manifold  of  the  nth  principal  quantum  number  due  to  the  presence 
of  an  ion  microfield  (ion  Stark  splitting)  is  given  by 

L$(e)D(u,  e)  = ^[d  • e,  D( u>,  e)] . (2.25) 


The  operator  Mj(e,  e']u)  describes  the  effects  of  the  motion  of  the  perturbing  ions  on 
the  radiator  dipole  through  fluctuations  in  the  ion  microfield  [38].  Using  an  ion  motion 
assumption  in  which  the  transitions  among  different  field  values  are  instantaneous  and 
uncorrelated  (Poisson  process),  the  ion-dynamics  collision  term  can  be  written  in  the 
form 

J de  Mi(e,  e ',u>)D(u>,  e ) = i upe>r  J de'Q(e')  ^p~^D(u,  e)  - p~}rD(u,  e')} 

= ivD(uj,e)  - iup€>r  f d e' Q(e')p~} rD{u,e)  (2.26) 


where  u is  directly  proportional  to  the  momentum  autocorrelation  function  of  the  ions 
and  is  a measure  of  the  importance  of  ion  motion  effects  on  the  line  shape.  Utilizing 
this  form  of  the  ion-dynamics  collision  term,  the  operator  D(u,  e)  can  be  shown  to  have 
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the  form 


where 


D(u,e ) = G(u,e) 


G(u>,  e)  — 


1 + w f de'  Q(e')G(u> , e’) 


/\  Pe,r  d 


u-  Lr-  L>J(e)  - Be(e)  - Me(u>,  e)  - \v 
and  the  distribution  function  has  the  form 


(2.27) 

(2.28) 


e-[Hr+d  e]/kBT 
Pe’r  = Trr{e-[^+d'dABT}  • 


(2.29) 


Finally,  we  arrive  at  the  main  topic  of  this  dissertation:  the  electron-broadening 
effects  of  the  perturbing  plasma  electrons.  The  effects  of  plasma  electron  perturbations 
on  the  radiator  dipole  moment  are  contained  in  the  static  part  Be(e)  and  the  dynamic 
part  Me(u,e)  of  the  electron-broadening,  or  width-and-shift  operator.  The  approxi- 
mation that  there  are  no  dynamical  correlations  between  the  electrons  and  the  ions  is 
reasonable  because  the  ions  are  much  heavier  and  slower  than  the  electrons  and  thus 
interact  with  the  radiator  dipole  over  a time  scale  much  longer  than  the  time  scale  of 
the  interaction  of  the  electrons  with  the  radiator  dipole.  The  electrons  affect  the  ion 
broadening  only  through  the  screened  ion  field  Esi.  The  width-and-shift  operator,  when 
expanded  out  in  a series,  has  the  form  of  a Born  expansion  in  the  electron-radiator 
interaction  Ve^r  ■ With  the  substitution  V$  — * A Ve[l\  the  expansion  has  the  form 


e)  = ]T  An{£(n)(<0  + e)}  (2.30) 

n = 1 


where  n defines  the  term  to  nth  order  in  the  interaction.  To  see  the  form  of  the  general 
width-and-shift  operator,  look  at  Equations  B.l  and  B.2  in  Appendix  B. 

In  order  to  work  with  a numerically  tractable  form  of  the  dynamic  part  of  the  width- 
and-shift  operator,  M’(w,  e),  a few  approximations  must  be  made.  These  approxima- 
tions are  discussed  in  detail  in  Appendix  B.  The  first  is  the  one-electron  approximation, 
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where  only  one  strong  electron  collision  is  allowed  with  the  radiator  at  a time.  As  a 
result  of  this  approximation,  correlations  between  plasma  electrons  have  been  neglected. 
Therefore,  the  screening  effects  of  electron  correlations  on  the  electron-radiator  inter- 
action are  approximated  with  a static  Debye  screening  model  [39],  or  in  other  words, 
interactions  between  the  radiator  and  the  perturbing  electron  are  considered  only  when 
the  electron  enters  a sphere  with  a radius  of  the  Debye  length  centered  around  the 
radiator. 

Next,  we  apply  the  no-quenching  approximation  where  no  perturbation- induced  non- 
radiative  transitions  are  allowed  between  upper  and  lower  states.  From  Equation  2.27, 
it  is  observed  that  the  width-and-shift  operator  acts  on  the  radiator  dipole  moment 
d.  Under  the  no-quenching  approximation,  only  matrix  elements  between  initial  and 
final  states  of  the  radiator  transition  for  this  occurrence  of  the  radiator  dipole  moment, 
noted  by  dj/,  are  allowed  in  the  line  shape  expression,  and  the  matrix  elements  of  the 
electron-radiator  interaction  between  initial  and  final  states  of  the  transition  are  zero. 
This  approximation  leads  to  four  separate  terms  in  the  width-and-shift  operator  as  is 
seen  in  Equation  B.17.  The  first  term  describes  the  broadening  on  the  upper  states  of 
the  transition,  the  last  term  describes  the  broadening  on  the  lower  states  of  the  tran- 
sition, and  the  two  intermediate  terms  are  interference  terms  between  upper  and  lower 
states.  Because  the  results  shown  in  this  work  consist  of  line  spectra  which  have  no 
significant  lower-state  broadening,  we  apply  a no  lower-state  broadening  approximation 
where  we  concentrate  specifically  on  the  first  term  that  includes  the  effects  of  upper-state 
broadening. 

With  these  approximations,  the  first  term  of  the  static  part  of  the  width-and-shift 
operator  has  the  form 


Bfirst(e)dif  = neTrle{V1 


le,r  P\e,e,r}  Pe,r^if 


(2.31) 
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and  the  first  term  of  the  dynamic  part  of  the  width-and-shift  operator  has  the  form 

r °o  - 

Mrir“(w,e)dl,  = -in,  J dte‘^H-'*“),hu{vffrlt)U(t,0) 

■ Pu,.,r  P,:i  d„  U<(t, 0)}  <?H-‘ . (2.32) 

The  time  development  operator  has  the  form 

U(t,0)  = ^exp|-i^tdt'yi{e1))r(f')|  (2.33) 

where  i)  indicates  a time-ordered  calculation  of  the  exponential  term, 

V$r(t)  = ei{H^+H^]t  Vffr  e~i{H^  r+H^1  (2.34) 

is  the  perturbing  electron-radiator  interaction  in  the  interaction  representation,  and 

H£tr  = Hr  + de.  (2.35) 

The  time  development  operator  satisfies  the  differential  equation 

= (2.36) 

Because  of  the  no-quenching  approximation,  operators  to  the  left  of  d^f  act  on  the  initial- 
state  manifold  of  the  transition  and  operators  to  the  right  of  d*/  act  on  the  final-state 
manifold  of  the  transition. 

The  form  of  the  Hamiltonian  that  now  defines  the  line  broadening  problem  as  pre- 
sented here  is  given  by 

H = H^T  + H[e  + V^T. 


(2.37) 
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Note  that  the  monopole  term  of  the  ion-radiator  interaction,  screened  by  the  plasma 
electrons,  has  been  incorporated  into  the  calculation  of  the  ion  electric  microfield.  The 
remainder  of  the  ion-radiator  interaction  is  approximated  solely  by  the  dipole  term  of 
the  interaction  that  depends  on  the  internal  state  of  the  radiator  and  the  ion  microfield. 
The  plasma  electrons,  in  the  presence  of  a charged  radiator,  are  described  by  the  plasma 
electron  Hamiltonian  and  the  monopole  term  of  the  interaction  between  the  perturbing 
electron  and  radiator  charge 

K = Ku  + V,'0’  (2.38) 

and  the  effects  of  electron  correlations  Ve<e  are  approximated  using  the  static  Debye 
screening  approach  mentioned  earlier. 

In  summary,  the  model  constructed  to  describe  plasma  broadened  line  shapes  include 
these  broadening  mechanisms: 

• Doppler  broadening  due  to  thermal  motion  of  the  radiator  is  included. 

• Stark  broadening  due  to  the  distribution  of  static-ion  fields  is  included. 

• Dynamical  ion  broadening  due  to  transitions  between  the  static-ion  fields  is 
included. 

• Electron  broadening  due  to  strong  electron-radiator  collisions  is  included. 

In  the  next  chapter,  a detailed  discussion  of  the  width-and-shift  operator  including 
the  effects  of  strong  collisions  is  given.  A classical  path  treatment  is  imposed  on  the 
perturbing  electrons,  and  a full-Coulomb  treatment  of  the  electron-radiator  interaction 
is  employed. 


CHAPTER  3 

ELECTRON-BROADENING  MODEL 


In  this  chapter,  the  electron  width-and-shift  operator  M'(uj,e)  is  discussed  in  more 
detail.  We  begin  with  a discussion  of  the  electron-radiator  interaction  used  in  the  calcu- 
lation of  M'(u>,e).  A classical-path  approximation  is  applied  to  the  perturbing  plasma 
electrons,  and  the  final  calculational  form  of  the  width-and-shift  operator  is  presented. 


3.1  Electron- Radiator  Interaction 

The  Hamiltonian  associated  with  line  broadening  as  shown  in  Chapter  2 has  the  form 

H = Hr  + H'le  + d-  e + Vffr  (3.1) 


where  Hr  is  the  isolated  radiator  Hamiltonian,  H[e  is  the  plasma  electron  Hamiltonian 
in  the  presence  of  the  monopole  term  of  the  electron-radiator  interaction,  d ■ e defines 
the  Stark  splitting  of  the  radiator  levels,  and  V^r  gives  the  interaction  of  the  perturbing 
electron  with  the  internal  electronic  structure  of  the  radiator.  The  monopole  term  of 
the  electron-radiator  interaction  has  the  form 

(0)  (Z-iV.)e2 

'le.r 


and  the  term  including  the  interaction  of  the  perturbing  electron  with  the  internal  elec- 
tronic structure  of  the  radiator  along  with  the  higher-order  multipole  terms  has  the 
form 


N, 


1 = 1 


( 

\>le  - ri 


(3.3) 
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where  Z is  the  atomic  number  of  the  radiator,  Nr  is  the  number  of  bound  electrons  on 
the  radiator,  and  xle  and  xl  are  the  positions  of  the  perturbing  electron  and  the  ith 
radiator  electron,  respectively,  in  the  radiator’s  center  of  mass  reference  frame. 

Under  a classical  treatment  of  the  perturbing  electron,  the  plasma  electron  Hamil- 
tonian including  the  monopole  interaction 


Hi  = Ku  + (3.4) 

describes  a perturbing  electron  as  a classical  particle  moving  in  a hyperbolic  trajectory 
as  it  passes  by  the  radiator.  Under  a quantum  mechanical  treatment  of  the  perturbing 
electron,  the  eigenvalue  problem  for  the  plasma  electron  Hamiltonian  is  defined  by 

h2k2 

Hle\k,l,m)  = — 1 k,l,m)  (3.5) 

zme 

where  | k.  I , m)  represents  a partial  wave  expansion  of  the  plasma  electron  state,  hk  is  the 
momentum  of  the  plasma  electron  state,  and  the  l and  m quantum  numbers  represent 
the  angular  momentum  of  this  state.  In  the  perturbing  electron  coordinate  representa- 
tion, the  states  ( x\k,l,m ) are  represented  by  Coulomb  wavefunctions.  The  differences 
between  the  classical  and  quantum  mechanical  treatments  are  discussed  further  in  the 
next  section. 

Now,  let  us  discuss  the  calculation  of  the  interaction  of  the  perturbing  electron  with 
the  internal  electronic  structure  of  the  radiator.  This  interaction  is  calculated  using  an 
expansion  in  terms  of  reduced  tensor  operators 

Nr  oo  t 
i — 1 t = 0 q = — t 


(3.6) 
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where 

A‘( r„r„)  = ^ and  C«(M)  = (j^l)  WM)  ' (3'7> 

The  t=0  term  represents  a monopole  contribution  between  the  perturbing  electron  and 
the  radiator  electron  for  perturbing  electrons  which  penetrate  the  orbitals.  The  largest 
contributions  to  this  interaction  come  from  the  t = 0 and  the  dipole  (t  = 1)  terms  of  this 
expansion  with  relatively  small  contributions  coming  from  the  higher-order  terms.  Due 
to  angular  momentum  selection  rules,  only  a finite  number  of  terms  in  this  expansion 
give  a nonzero  contribution  to  the  interaction.  Even  though  the  monopole  term  has  been 
removed  from  this  expansion,  a calculation  including  all  of  the  terms  in  the  expansion  is 
referred  to  throughout  the  remainder  of  this  dissertation  as  a full-Coulomb  calculation. 
A calculation  including  only  the  t = 1 term  is  a part  of  the  dipole  approximation,  and 
the  dipole  approximation  is  discussed  in  more  detail  in  the  first  section  of  Chapter  4. 

One  question  that  may  come  up  at  this  point  is  why  the  dipole  only  calculation  is  valid 
in  the  radiator-ion  interaction  but  not  valid  in  the  radiator-electron  interaction.  One 
must  remember  that  for  highly  charged  radiators,  the  perturbing  ions  experience  a strong 
repulsion  due  to  the  monopole  interaction  with  the  radiator  through  the  calculation  of 
the  ion  microfield,  and  thus  are  not  very  likely  to  penetrate  the  orbitals  of  the  radiator. 
Also,  the  terms  of  the  expansion  higher  than  the  dipole  term  have  been  shown  to  have 
a relatively  small  effect  [37]  in  comparison  to  the  electron-broadening  effects.  Of  course, 
this  is  definitely  not  the  case  for  the  perturbing  electron  which  experiences  a strong 
attraction  to  the  radiator  and  can  even  penetrate  the  radiator  orbitals,  thereby  resulting 
in  significant  contributions  from  the  t = 0 and  dipole  (t  = 1)  terms. 
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3.2  Calculation  of  the  Width-and-Shift  Operator 

As  was  shown  in  Chapter  2,  the  first  term  of  the  static  part  of  the  width-and-shift 
operator  has  the  form 


Bfirst(e)dif  = neTrle{V^rplete,r}p-ldif 


(3.8) 


and  the  first  term  of  the  dynamic  part  of  the  width-and-shift  operator  has  the  form 


M 


r oo 

fiTat{u ,,e)dif  = -i  neJ  dte^-H^+i^Tvu[vffr(t)U(t,0) 


■ Pie,  e, r VPr  P~l  dif  U^t,  0)}  ^ . (3.9) 


The  time  development  operator  has  the  form 


U {t,  0)  = exp  | -i  d t'  Vffr 


(0 


(3.10) 


where  d indicates  a time-ordered  calculation  of  the  exponential  term, 


V^r(t)  = eH We, r +//(,)<  ydO  e~'{Hl’r+H^]t 


(3.11) 


is  the  perturbing  electron-radiator  interaction  in  the  interaction  representation,  and 


He  r — Hr  -f-  d • 6 . 


(3.12) 


The  time  development  operator  satisfies  the  differential  equation 


(313) 


Note  that  for  large  t , the  matrix  elements  of  U (t,  0)  become  the  5'-matrix  for  electron 
scattering  by  the  radiator  [23].  In  the  context  of  this  theory,  these  completed  collisions 
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give  the  impact  limit  for  small  frequency  separations  near  line  center.  For  short  times, 
U(t,  0)  approaches  1,  and  an  expansion  of  the  radiator  resolvent,  represented  by  equa- 
tion 2.28,  in  terms  of  the  width-and-shift  operator,  Be(e)  -f  Me(u>,  e),  leads  to  an  appro- 
priate form  for  the  high  frequency  wings  of  the  line  shape  [24],  The  theory  presented 
here  is  therefore  “unified”  in  the  sense  that  it  describes  accurately  both  the  center  and 
the  wings  of  the  line  shape. 

Continuing  in  this  section,  a classical-path  approximation  is  applied  to  the  trace 
over  the  perturbing  electron  coordinates  in  the  width-and-shift  operator.  Because  the 
perturbing  electrons  move  in  hyperbolic  trajectories,  a rotational  transformation  is  also 
performed  from  an  atomic  collision  reference  frame  to  a hyperbolic  collision  reference 
frame  in  order  to  simplify  the  calculation  of  the  matrix  elements  of  the  width-and-shift 
operator.  We  end  this  section  by  writing  out  the  final  expression  for  the  calculation 
of  the  matrix  elements  of  the  width-and-shift  operator  and  the  matrix  elements  of  the 
interaction 

3.2.1  Classical-Path  Approximation 

Up  to  this  point,  the  treatment  of  the  perturbing  electrons  as  either  quantum  mechan- 
ical or  classical  particles  has  not  been  specified.  In  this  section,  an  argument  is  made 
that  a classical  treatment  of  the  perturbing  electrons  is  appropriate  in  the  calculation 
of  the  width-and-shift  operator  if  an  appropriate  length  scale  is  applied  to  the  problem. 

In  applying  the  classical-path  approximation,  the  wave  packets  of  the  perturbing 
electrons  that  are  not  interacting  with  the  radiator  are  assumed  to  be  small  enough  so 
that  they  do  not  significantly  overlap  with  each  other  [21].  Also,  when  the  perturbing 
electrons  interact  with  a highly  charged  radiator  such  as  Ar+16  or  Ar+17,  there  is  a large 
acceleration  of  the  electron  as  it  passes  by  the  radiator,  thus  resulting  in  a smaller  wave 
packet  size  in  the  region  of  the  interaction.  The  decrease  in  the  size  of  the  wave  packet 
due  to  the  acceleration  and  the  very  short  time  in  which  the  electrons  remain  in  close 
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proximity  to  the  radiator  allow  for  a classical  treatment  of  the  perturbing  electrons, 
even  those  which  penetrate  the  orbitals  of  the  radiator,  as  long  as  we  apply  a correction 
due  to  quantum  diffraction  effects.  This  correction  is  incorporated  through  a minimum 
distance  cutoff  which  is  discussed  later  in  this  section  and  also  in  Appendix  C.  One 
possible  drawback  to  the  smaller  wave  packets  is  that  they  do  not  hold  together  very 
long,  but  the  electrons  are  also  moving  more  quickly  near  the  radiator  and  thus  should 
hold  together  at  least  during  the  time  of  collision. 

It  has  been  previously  shown  that  the  classical-path  approximation  is  appropriate 
for  both  the  neutral  and  charged  cases  in  the  VCS  model  when  [21] 


V, 


(i) 


le,  r 


ih 


< l 


(3.14) 


where  r is  the  duration  of  the  collision  and  lh  is  the  angular  momentum  of  the  perturbing 
electron.  This  is  also  the  condition  that  defines  a strong  collision  in  their  model  if  l is  on 
the  order  of  one  or  smaller.  However,  this  expression  is  appropriate  only  for  calculations 
in  which  the  monopole  term  due  to  the  penetration  of  the  orbitals  is  not  included  in  the 
calculation. 

For  the  case  of  the  highly  charged  radiator  where  the  monopole  term  due  to  the 
penetration  is  included,  the  interaction  strength  changes  dramatically  over  the  region  of 
the  interaction.  However,  this  is  offset  by  the  significant  acceleration  of  the  electron  near 
the  charged  radiator  as  the  electron  moves  in  a hyperbolic  path  because  the  electron 
spends  less  time  in  the  region  of  the  strongest  interaction.  Because  the  speed  of  the 
electron  is  affected  by  the  charge  of  the  radiator  Z and  the  temperature  of  the  electrons 


in  the  plasma  k gT,  the  largest  value  of  the  interaction 


Vi 


(i) 


1 e,r 


that  we  allow  should 


be  related  to  these  quantities.  The  typical  temperature-  and  Z- dependent  minimum- 
distance  cutoff  applied  to  the  minimum  separation  of  the  perturbing  electron-charged 
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radiator  collision  is  rmin  = Z riand , where  riand  is  the  Landau  length 


T land 


(3.15) 


This  type  of  cutoff  sets  the  maximum  of  of  the  electron-radiator  interaction  to  approxi- 
mately 


(i) 


kRT 


(3.16) 


and  assures  that  the  strength  of  the  interaction  is  balanced  by  the  speed  of  the  electron 
as  it  passes  by  the  charged  radiator.  In  Appendix  C,  we  derive  a cutoff  of  the  form 


Z V land 

rmin  = In  9lJ(r  ->  0)  • 

Note  that  \ngij(r  — > 0)  modifies  the  cutoff  for  quantum  diffraction  effects  between  the 
perturbing  electron  and  the  radiator  [40]  because  there  is  overlap  between  the  wave 
packet  of  the  perturbing  electron  and  the  radiator  orbitals.  With  the  limit  on  the 
strength  of  the  interaction  and  the  inclusion  of  quantum  corrections  through  the  cutoff, 
the  classical-path  approximation  can  be  used  for  perturbing  electrons  interacting  with 
highly  charged  radiators. 

Under  a classical-path  approximation,  there  are  three  important  changes  made  in 
the  expression  for  the  width-and-shift  operator.  The  first  change  involves  replacing 
the  portion  of  the  density  matrix  associated  with  the  perturbing  electron  Hamiltonian 
with  its  classical  analogue.  This  includes  the  treatment  of  the  velocity  distribution  of 
the  plasma  electrons  outside  the  region  of  interaction  as  being  in  thermal  equilibrium 
through  a Maxwell-Boltzmann  distribution.  This  also  implies  that  electron  degeneracy 
effects  are  small,  thus  requiring  the  condition 


E Fermi  kgT 


(3.18) 
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where 

E Fermi  = {37 r2ne}2/3  Rydbergs  (3.19) 

is  the  Fermi  energy  for  a plasma  of  density  ne. 

The  second  change  involves  the  replacement  of  the  quantum  mechanical  expression 
for  the  time  dependent  electron-radiator  interaction  with  an  expression  for  the  interac- 
tion which  depends  on  the  classical  time-dependent  position  and  velocity  of  the  perturb- 
ing electron  in  its  hyperbolic  trajectory  about  the  radiator.  In  terms  of  an  expression, 
this  replacement  has  the  form 

j dx  J dx  elh^k  2~k^t/2me  (k'  ,l'  ,m'\x')(x'\V^r(rad)\x)(x\k,l,rn) 

— ► V$r(rad,xle(t),vie(t))  (3.20) 

where  “rad”  indicates  the  dependence  on  the  internal  states  of  the  radiator  and  | k , l,  m) 
represents  a partial  wave  expansion  of  the  perturbing-electron  state.  As  stated  in  the 
previous  section,  these  states  are  the  appropriate  eigenfunctions  for  the  electron  Hamil- 
tonian H[e 

h2k2 

Hu\k,l,m)  = -^-\k,l,m) , (3.21) 

and  in  the  perturbing  electron  coordinate  representation,  the  states  (a;|  k,  l , m)  are  rep- 
resented by  Coulomb  wavefunctions. 

The  third  change  involves  the  calculation  of  the  one-electron  trace  over  the  free 
electron  perturber  states.  In  the  quantum  mechanical  calculation,  the  one-electron  trace 
is  written  as  an  integral  over  momenta  states  and  a sum  over  angular  momenta  states 

°°'  1 r kDeb 

Trie{- • •}  = 2^  2_2  / dk  k2(k  ,m  ,l\ • ■ ■ \k  ,1  ,m) . 

1=0  m=—l 


(3.22) 
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Under  the  classical-path  approximation,  the  one-electron  trace  for  the  electron  trajectory 
has  the  form 


The  variables  (speed  outside  interaction  region),  b (impact  parameter),  and  Q (ori- 
entation of  trajectory)  define  an  unique  hyperbolic  path  about  the  radiator.  The  lower 
cutoff  on  the  impact  parameter  is  set  so  that  the  minimum  distance  of  approach  between 
the  electron  and  the  radiator  is  given  by  Equation  3.17.  In  order  to  account  for  the 
removal  of  electron  correlations,  the  upper  limit  on  the  impact  parameter  is  set  so  that 
electrons  outside  a sphere  of  radius  the  Debye  length  are  effectively  screened  from  the 
radiator.  In  order  to  define  the  initial  position  of  the  electron  on  the  hyperbola,  we 
also  need  to  integrate  over  the  time  to  that  it  takes  the  electron  to  leave  the  interaction 
region  while  moving  from  any  point  within  this  region  along  its  hyperbolic  trajectory. 
In  other  words,  we  sum  up  the  contributions  from  all  possible  incomplete  collisions.  The 
upper  limit  tcoi  is  the  total  time  that  it  takes  an  electron  to  move  through  the  interaction 
region  along  its  trajectory.  Note  that  during  the  time  df0,  the  number  of  electrons  that 
begin  a collision  with  the  radiator  is  given  by  ne(27rbdb)voodt0-  Also,  the  integration 
over  this  electron  number,  multiplied  by  the  velocity  distribution  of  the  plasma  electrons 
outside  the  interaction  region,  gives  the  appropriate  average  over  all  possible  collisions 
of  the  free  electrons  with  the  radiator  for  a given  orientation  of  the  electron  trajectories. 

Throughout  the  rest  of  this  work,  the  classical  description  of  the  perturbing  electrons 
is  used,  including  the  classical  forms  of  the  electron-radiator  interaction  and  the  one- 
electron  trace. 

3.2.2  Radiator  Basis  and  the  Matrix  Elements  of 


A complete  set  of  basis  states  must  now  be  chosen  for  the  isolated  radiator  in  order 
to  calculate  the  matrix  elements  of  the  width-and-shift  operator.  Because  the  radiator 
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species  considered  in  this  dissertation  must  be  treated  relativistically,  the  natural  energy 
basis  for  the  isolated  radiator  is  the  total  angular  momentum  energy  basis  | a , j , m ) , 
where  a represents  all  of  the  quantum  numbers  excluding  j and  m.  Primes  on  the 
quantum  numbers  refer  to  a lower  or  final  state  of  the  transition. 

Because  the  classical  trajectories  of  the  perturbing  electrons  in  the  charged-radiator 
problem  are  hyperbolic,  it  is  useful  to  perform  a rotational  transformation  from  an 
atomic  collision  reference  frame  in  which  the  microfield  e points  along  the  z-axis  to  a 
hyperbolic  collision  reference  frame  when  calculating  the  matrix  elements  of  the  width- 
and-shift  operator.  In  the  hyperbolic  frame,  the  x-axis  is  aligned  so  that  it  points  outward 
from  the  radiator  toward  the  point  of  closest  approach  of  the  perturbing  electron  in  its 
hyperbolic  trajectory.  The  z-axis  points  outward  from  the  radiator  perpendicular  to  the 
plane  of  the  trajectory.  The  time  development  operator  in  this  hyperbolic  frame  is  given 


by  [25] 


Uh(t,0)  = D(Cl)U(t,0)D-\n) 


(3.24) 


where  D(Q)  is  the  rotation  operator  which  rotates  from  the  atomic  collision  frame  men- 
tioned above  to  the  hyperbolic  collision  frame.  This  time  development  operator  satisfies 
the  differential  equation 


ihdUh^°)  =vh(t)uh(t,o) 


(3.25) 


where 


v»(i)  = D(n)V™(t)D-'((l) 


(3.26) 


is  the  electron-radiator  interaction  as  seen  in  the  hyperbolic  frame. 


Using  the  total  angular  momentum  basis  set  and  the  rotational  transformation  men- 
tioned above,  the  matrix  elements  of  the  static  part  of  the  width- and-shift  operator  have 
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the  form 


r oo 

{<x'aj'a,m'a'’aa,ja,rna\Bfirst(e)\a'b,j'b,m'b]ab,jb,mb)  = ne  / dv^An 

J 0 

fbr>eb  rtc°i  f df2  . 

X / db2nb  / dto^oc  ^ J {^m00,m!(«a,ia,mc  | 

J bmin  ^ 0 J mm. 


m c,  m d 


x^(0)p!t,rpS 


(3.27) 


and  the  matrix  elements  of  the  dynamic  part  of  the  width-and-shift  operator  have  the 
form 


n oo 

(a'aJa,m’a',aaJa,rna\Mfirst(uj,e)\a'b,j'b,m'b-,ab,jb,mb)  = -me  / dt 

Jo 

/OO  1-bDeb  rtcol  f HO  . 

d^oo47 rjf  d6  27T&  jf  dtoVoo  j — (4,™! 


77lc  j 771  d ; 777. c , 171^ 


x (aa,ia,"ic|e  lH^)rt/nVh(t)  Uh(t,0)  p[hJeir  Vh(0)  p£l  l\ab,jb,md) 
x DmimbD^~Ua'b,j'b  ,m'd  \Ul(t,  0)  e'^l  , j'a  ,m'c)D^?K  | (3.28) 


where 


P{eh)r  = D^p^Diny1 

(3.29) 

P^e,  r = DWpu^DW-1 

(3.30) 

H%1  = Hr  + D{tt){d  • e)D{ny1 

(3.31) 

Dm]m'  = (n’j’m\D(Q)\n’j’m')  ■ 

(3.32) 

The  only  terms  in  the  calculation  of  the  matrix  elements  of  the  static  and  dynamic 
width-and-shift  operator  which  are  dependent  on  the  orientation  of  the  trajectory  fl 
are  the  Stark-broadening  terms  d ■ e and  the  matrix  elements  of  the  rotation  operators 
D ^ ^ , . All  other  terms  are  spherically  symmetric. 

771  , 771  1 J u 
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To  simplify  this  expression  further,  an  approximation  is  needed  to  deal  with  the  SI 
dependence  of  the  Stark-broadening  term.  At  this  point,  we  remove  the  dependence  of 
the  width-and-shift  operator  on  the  ion  microfield  e,  as  was  done  in  other  models  [18], 
thereby  leaving  the  rotation  terms  ^ , outside  of  both  the  time  development  operator 
and  the  reduced  distribution  functions  as  the  only  terms  dependent  on  Q.  With  this 
approximation,  the  expressions 


/m  = f E)Ua)-i  £)Cn) 

\ /l—  I g^2  ma,mc  md,mb 


(3.33) 


nOo)-l  n(ii)  ndt)  1 nb'a) 
™a,mc  ^md,mb  ^ m'b,m'd  m’c,m'a 


(3.34) 


represent  a spherical  average  over  all  possible  orientations  of  the  electron  trajectories. 
The  calculational  details  of  these  averages  have  been  presented  in  previous  papers  by 
various  authors  [24,25].  Thus,  the  matrix  elements  of  the  width-and-shift  operator  under 
these  averages  have  the  form 

r oo 

(^'aj'a,rn'a]aa,ja,rna\Bfirat\a'b,j'b,m'b]ab,jb,mb)  = ne  / du^Tr 

Jo 

r &De6  r tool  1 % 

x / db2nb  d t0  -----  — Y,  { ( a a , 3 a , m c ] 14(0)  p Ut  r p jT1 

Jbmin  Jo  * Ja  + 1 mc 

x | <2 6 , j a imc)}  ^ja,jb^rna,rnbSa>aiafbSj^j^8mr^m^  (3.35) 


and 


r < 

(oi'a,ja,m'a]aa,ja,ma\Mfirst(uj)\a'b,j'b,rn'b]ab,jb,rnb)  = -ine 

Jo 
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(3.36) 
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If  final-state  broadening  effects  are  to  be  left  out  of  the  calculation,  the  first  term  of 
the  dynamic  part  of  the  width-and-shift  operator  must  be  modified  to  maintain  consis- 
tency by  replacing  U\t,  0)  with  the  unity  operator.  Therefore, 


E E (3-37) 

and  the  dynamic  part  simplifies  in  a manner  similar  to  that  of  the  static  part.  The  first 
term  of  the  static  part  does  not  contain  any  final-state  broadening  contributions  and 
remains  unchanged.  Therefore,  the  matrix  elements  of  the  static  and  dynamic  parts  of 
the  width-and-shift  operator  are  now  given  by 


(<* a Ja  ,ma-,aa  ,j  a ,ma  \ Bimtial\ab  ,jb  ,mb;ab  ,jb  ,mb ) 
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dwoo  47 r 


r^Deb  rtcol  2 ^ 

x/  db2nb  d to  v ^ - -7-  -■  V"  { ( aa  , ja  , mc  \ Vh{0)  ple,r  P r* 

Jbmin  Jo  Ja.  m ^ 

X | ab  > ja  ) }}  ^ja,jb^rna,mb^a'aia'b^j'aj'b^m'aim'b  (3.38) 


and 


r oo 

(^aJ'a^m'a-,aa,ja,ma\Mfirst(uj)\a'bj'b,m'b-,ab,jb,mb)  = -ine  / d t 

Jo 

rf>Deb  rtcol  1 
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Jo  Jbmin  Jo  ZJa  + 1 mc 

x e-i  Hrt/hyh^  Uh^  pier  Vfc(0)  p~x\  ab , ja  ,mc)(a'a,j'a,m'a  \ eiHrt/n 
X I aa  ij a )ma)}  $ja,  jfc^m0,m6^a'a>Q'  ^ j^m'a,m'b  ' (3.39) 


To  calculate  the  matrix  elements  of  the  width-and-shift  operator  for  initial-state 
broadening,  matrix  elements  of  the  time  development  operator  and  of  the  electron- 
radiator  interaction  Vh(t)  over  the  initial  states  are  needed.  To  preserve  the  appropriate 
time  ordering  of  the  interaction  terms  in  the  expansion  of  the  time  development  operator, 
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a coupled  set  of  first-order  differential  equations  of  the  form 

d 

ih—{aa,ja,ma\Uh{t,0)\ab,jb,mb)=  ^ {aa  ,ja  ,ma\Vh(t) 

etc  Jc,  rnc 

x\ac,jc,mc)(ac,jc,mc\Uh{t,0)\ab,jb,mb)  (3.40) 

must  be  solved.  The  number  of  differential  equations  in  this  coupled  set  is  TV2,  where 
N is  the  number  of  radiator  states  included  in  the  initial-state  manifold.  Many  of  these 
matrix  elements  are  zero  because  of  angular  momentum  selection  rules  in  the  calculation 
of  the  matrix  elements  of  Vb(t). 

As  a matter  of  convenience,  another  reference  frame  can  be  defined  in  which  the 
x-axis  moves  with  the  electron  as  it  traverses  its  hyperbolic  trajectory  past  the  radiator. 
The  rotation  operator  that  defines  the  rotation  from  the  hyperbolic  reference  frame  to 
this  rotating  reference  frame  is  be  given  by 

D{(f>ie{t))  = e-^u{t)Jz/h  (3.41) 

where  4>ie(t)  is  the  azimuthal  angle  and  Jz  is  the  zth  component  of  the  total  angular 
momentum  operator.  With  this  rotation,  the  matrix  elements  of  the  electron-radiator 
interaction  in  the  hyperbolic  frame  can  be  written  in  the  form 

( ua,ja,ma  \Vh(t)\ab,jb,mb)  = e^a-m^ieW 

x(<*a,  ja  ,ma  \Vrot(t)\ab,jb,mb)  (3.42) 

where  Vrot(t)  is  the  electron-radiator  interaction  in  this  rotating  reference  frame.  Keeping 
in  mind  that  the  electron  coordinates  are  being  treated  classically,  the  radiator  matrix 
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elements  of  this  interaction  can  be  given  by 


OO 


{oia,ja,ma\Vrot(t)\ab,jb,mb)  = {-l)u  m“ 


ja  t jb 


(3.43) 


i = 1 


\-ma  q mb  j 


with  the  application  of  the  Wigner-Eckart  theorem  [41,42].  Note  that  in  this  rotating 
reference  frame,  9\e{t)  = 7t/2  and  (j)ie{t)  = 0. 

The  reason  for  writing  the  potential  in  terms  of  irreducible  tensor  operators  and  then 
applying  the  Wigner-Eckart  theorem  is  now  clear.  The  selection  rules  for  the  occurrence 
of  nonzero  matrix  elements  of  the  interaction  are  now  embodied  in  the  calculation  of  the 
3J  symbol 


ja  t jb 


(3.44) 


\-ma  q mb  I 

which  is  zero  unless  the  triangle  relations 


t ^ I ja-jb 


(3.45) 


t ^ ja+jb 


(3.46) 

(3.47) 

(3.48) 


q = ma  — mb 
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hold.  Further  properties  of  the  3J  symbols  are  discussed  in  Appendix  D.  Therefore, 


ja+jb 

{oia,ja,ma\Vrot{t)\ab,jblmb}  = ^ 

t =max[\ma-mb\,\ja-jb\] 
Nr 


X 


t = l 

ja  t 

-ma  ma  — mb  mb 


. \ 

Jb 


(3.49) 


J 


The  line  shapes  to  be  presented  in  this  dissertation  result  from  transitions  in  hydro- 
genic  radiators.  Therefore,  the  next  step  is  to  uncouple  the  spin  angular  momenta  from 
the  orbital  angular  momenta.  Because  the  interaction  Vrot(t ) does  not  depend  on  the 
spin  coordinates,  the  matrix  elements  of  this  interaction  have  the  form 


(aa,ja,ma \Vrot(t)\ab,jb  ,mb)  = 6SatSb(-  1)J  ° m°+:,6+So  ( [j  a,jb\[la,lb\)1/2 

min  [ja+jb,la+h] 


t =max[  \ma  — Tnb\,\j a—j lb\] 
j a t 


/ 

la  t lb 

o 

o 

V 

—ma  ma-mb  mb 


J b 1 . | l a s a j a 

(-1)  . 

j b t lb 


where  the  last  term  is  a 6J  symbol  and  the  definitions 

[laJb,"']  = (21  a + l)(2lb  + 1) 


(3.50) 


(3.51) 


(na,/a||At(rie)C'i(t)||n6,/6)  = (-1)‘“[Z„,  h}1'2 


l a t 


0 0 0 
x(na,/a||A  (ru)\\nb,lb) 


V 


(3.52) 
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have  been  used.  The  definition  of  the  6J  symbol  is  given  in  Appendix  D,  along  with 
the  specific  treatment  of  the  6J  symbol  occurring  in  this  expression.  Also,  note  that  the 
uncoupling  of  the  orbital  and  spin  angular  momenta  places  a third  constraint  on  the 
summation  over  t.  Finally,  the  reduced  matrix  element  in  the  tth  term  of  the  multipole 
expansion  is  given  by 


(na,/a||At(rie)||n6,/6) 


dr  r2  R] 


» 


,t+ 1 


$t,  o 

rie 


Rnb,ib(r)  (3-53) 


where  Rnb,lb(r ) are  the  relativistic  radial  wave  functions. 

In  the  next  chapter,  we  begin  by  discussing  the  computational  methods  used  in 
calculating  the  width-and-shift  operator  described  in  the  current  chapter  and  its  incor- 
poration into  the  calculation  of  theoretical  spectral  lines.  We  then  discuss  the  various 
effects  of  the  width-and-shift  operator  on  the  shape  and  location  of  the  theoretical  spec- 
tral lines.  Comparisons  are  then  made  between  calculated  spectral  lines  utilizing  the 
electron-broadening  model  described  in  this  chapter  and  calculated  spectral  lines  utiliz- 
ing other  electron-broadening  models.  The  process  of  line  merging  is  also  discussed. 


CHAPTER  4 

COMPUTATIONAL  ANALYSIS 

In  this  chapter,  the  computational  methods  and  results  of  the  semiclassical  (classical- 
path)  model  to  all-order  in  the  full-Coulomb  electron-radiator  interaction  as  described 
in  the  previous  chapters  are  discussed  in  detail.  The  dependencies  of  line  widths  and 
shifts  on  various  plasma  and  radiator  state  parameters  are  considered.  Comparisons  of 
calculated  line  shapes  using  full-Coulomb  and  dipole-approximation  versions  of  the  semi- 
classical all-order  model  are  shown  to  demonstrate  the  importance  of  the  full-Coulomb 
treatment  of  the  electron-radiator  interaction.  We  then  compare  calculated  lines  shapes 
using  full-Coulomb  and  dipole-approximation  versions  of  both  the  semiclassical  all-order 
model  and  the  quantum  mechanical  model  to  second-order  in  the  electron-radiator  inter- 
action. We  show  that  the  full-Coulomb  all-order  and  second-order  models  in  the  electron- 
radiator  interaction  differ  mainly  in  the  wings  of  the  line  shape  and  in  the  calculation 
of  merging  lines.  Note  that  all  line  shapes  shown  throughout  this  chapter  are  calculated 
with  an  atomic  fraction  concentration  of  1%  Argon  in  D 2 because  the  width- and-shift 
operator  does  not  depend  on  this  concentration. 

4.1  Computational  Methods 

The  computation  of  spectra  has  been  broken  into  two  parts.  The  first  part  involves 
the  calculation  of  matrix  elements  of  a width-and-shift  operator  using  a code  based 
on  the  expressions  presented  in  the  previous  chapter.  The  width-and-shift  values  are 
incorporated  into  a second  code  which  uses  these  values  to  generate  theoretical  spectra. 
The  theoretical  spectra  can  then  be  used  to  analyze  experimental  spectral  data.  This 
section  concentrates  on  the  computation  of  line  widths  and  shifts  and  the  corresponding 
theoretical  spectra. 
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To  calculate  the  line  widths  and  shifts  and  the  line  profiles,  information  concerning 
the  the  relativistic  atomic  physics  must  first  be  incorporated  into  the  computational 
codes.  The  line  shape  code  and  the  width-and-shift  code  written  by  the  author  of  this 
dissertation  require  as  input  the  relativistic  energy  level  structure  generated  by  a suite  of 
atomic  physics  codes  written  by  Robert  Cowan  at  Los  Alamos  National  Laboratory  [43]. 
In  addition,  relativistic  radial  wavefunctions  for  the  bound-electron  states  are  extracted 
from  this  suite  of  atomic  physics  codes  to  be  used  in  the  width-and-shift  code. 

The  calculation  of  the  width-and-shift  values  begins  with  the  generation  of  the  matrix 
elements  of  the  full-Coulomb  electron-radiator  interaction  in  the  rotating  reference  frame 
(aa  , ja  ,ma  |Kot(^)|  &b  , jb  ,mb)  as  described  in  the  previous  chapter.  The  advantage 
of  using  this  rotating  reference  frame  is  that  the  matrix  elements  are  not  dependent 
on  the  azimuthal  angle  (f)(t ) and  are  only  dependent  on  the  perturbing  electron’s  radial 
distance  from  the  center  of  the  radiator,  thereby  reducing  the  memory  needed  for  storage. 
Remember  that  in  this  rotating  reference  frame,  the  z-axis  is  perpendicular  to  the  plane 
of  the  trajectory,  and  the  x-axis  moves  with  the  perturbing  electron  as  it  traverses  its 
hyperbolic  trajectory.  The  purpose  of  calculating  the  matrix  elements  first  and  then 
storing  them  in  memory  is  that  it  is  much  faster  to  read  them  from  memory  then  it  is 
to  generate  them  as  they  are  needed  in  the  calculation  of  the  width-and-shift  values. 

To  calculate  these  matrix  elements,  the  reduced  matrix  elements  of  the  irreducible 
radial  tensor  operator  (na  ,la  ||^4f(rie)||  nb , h ) are  needed,  and  these  require  the  relativis- 
tic radial  wavefunctions  Rnj(r).  Although  the  radial  wavefunctions  for  the  hydrogenic 
radiator  problem  can  be  generated  easily,  the  radial  wavefunctions  used  in  this  model 
are  extracted  from  Cowan’s  suite  of  atomic  physics  codes.  This  approach  allows  for  an 
easier  expansion  of  this  code  to  handle  helium-like  radiators  where  the  calculation  of  the 
relativistic  radial  wavefunctions  are  no  longer  quite  as  simple. 

After  the  calculation  of  the  matrix  elements  of  the  electron-radiator  interaction,  the 
integrations  over  the  initial  electron  velocity  v^.  the  impact  parameter  6,  and  the  time 
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of  closest  approach  t0  are  carried  out  over  the  appropriate  limits  on  the  integrands  of 
both  the  static  B and  dynamic  M (u>)  parts  of  the  width-and-shift  operator.  In  addition 
to  the  dependence  of  these  integrands  on  the  matrix  elements  of  the  electron-radiator 
interaction,  the  integrand  of  the  dynamic  part  also  depends  on  matrix  elements  of  the 
time  development  operator.  The  placement  and  calculation  of  the  integrand  points 
are  weighted  more  heavily  in  those  regions  in  which  the  integrand  is  changing  most 
rapidly.  These  three  integrations  represent  the  one-electron  trace  as  defined  in  the 
previous  chapter. 

For  each  value  of  the  initial  electron  velocity  and  the  impact  parameter,  the  matrix 
elements  of  the  time  development  operator  need  to  be  calculated.  This  is  accomplished 
by  solving  a coupled  set  of  first-order  differential  equations  of  the  matrix  elements  of  the 
time  development  operator,  thereby  including  the  appropriate  time  ordering  in  the  time 
development  operator.  These  differential  equations  depend  specifically  on  the  behavior 
of  the  matrix  elements  of  the  electron-radiator  interaction.  Out  of  the  different  ordinary 
differential  equation  solvers  that  were  tested,  a 4th  order  Runge-Kutta  solver  adjusted  for 
complex  numbers  with  5th  order  error  checking  resulted  in  the  most  efficient  calculation 
of  this  set  of  first-order  differential  equations. 

The  final  step  in  the  calculation  of  the  width-and-shift  values  is  to  calculate  the 
Laplace  transform  in  the  dynamic  part  of  the  width-and-shift  operator.  This  transform  is 
relatively  simple  to  calculate.  Most  of  the  calculational  effort  in  obtaining  the  width-and- 
shift  values  occurs  in  the  calculation  of  the  one-electron  trace  in  the  dynamic  part  of  the 
width-and-shift  operator  due  to  its  dependence  on  the  time  development  operator.  The 
dynamic  part  generally  takes  about  100  times  longer  than  the  static  part  to  calculate. 

To  generate  the  theoretical  spectra,  a previous  line  shape  code  known  as 
MERL  [35,44]  is  used  with  the  exception  that  the  quantum  mechanical,  multi-electron, 
dipole-approximation  treatment  of  the  electron  width-and-shift  operator  is  replaced  with 
the  model  presented  in  this  dissertation.  In  order  to  function,  this  line  shape  code  needs 
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the  energy  level  structure  of  the  radiator.  It  also  needs  the  reduced  dipole  matrix  ele- 
ments that  are  used  in  the  calculation  of  the  radiator  dipole  autocorrelation  function 
and  the  Stark  splitting  term.  As  stated  above,  this  information  is  supplied  by  Cowan’s 
suite  of  atomic  physics  codes. 

If  one  wishes  to  calculate  the  width- and-shift  values  under  the  dipole  approximation, 
the  same  procedure  stated  above  for  calculating  the  width-and-shift  operator  can  be 
used  with  two  exceptions.  The  first  exception  is  to  keep  only  the  dipole  term  in  the 
multipole  expansion  of  the  electron-radiator  interaction  V^r-  The  monopole  contribu- 
tion Vie],  that  was  added  to  the  perturbing  electron  Hamiltonian  is  not  affected  by  this 
approximation.  The  second  exception  is  to  set  the  lower  limit  on  the  integration  over  the 
impact  parameter  to  zero  to  compensate  for  ignoring  the  monopole  orbital  penetration 
term.  Because  almost  all  of  the  static  shift  results  from  this  monopole  term,  there  is  no 
significant  shift  and  distortion  under  the  dipole  approximation.  This  classical  version  of 
the  dipole  approximation  is  analogous  to  the  quantum  mechanical  dipole  approximation 
used  in  MERL. 


4.2  Analysis  of  Line  Widths  and  Shifts 

In  this  section,  line  spectra  calculated  from  a semiclassical  model  to  all-order  in  the 
perturbing  electron-radiator  interaction  (SC  all-order  model)  is  presented  with  a focus  on 
spectra  from  hydrogenic  argon  radiators,  Ar+17,  immersed  in  the  plasma.  Dependencies 
of  the  width-and-shift  of  line  shapes  on  the  electron  density  and  temperature  and  on  the 
principal  and  orbital  angular  momentum  quantum  numbers  is  discussed  in  some  detail. 

In  general,  the  width  and  shift  of  spectral  lines  increases  with  increasing  principal 
quantum  number  n of  the  upper  states  of  a transition  in  a Rydberg  series.  The  physical 
reasoning  for  this  behavior  is  that  the  higher-lying  energy  shells  of  the  radiator  provide 
a larger  spatial  volume  for  the  interaction  of  the  perturbing  electron  with  the  corre- 
sponding radiator  state,  or  in  other  words,  the  higher  n states  are  less  tightly  bound 
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and  more  easily  perturbed.  Also,  the  sub-shells  with  lower  orbital  angular  momentum 
l provide  a larger  spatial  volume  for  the  interaction.  Thus,  the  spherical  l = 0 or  s 
orbital  experiences  the  largest  width  and  shift,  with  smaller  width- and-shift  values  as  / 
increases.  This  behavior  is  shown  in  Figure  4.1  for  width-and-shift  values  at  line  center 
with  an  electron  density  of  1 x 1024  cm-3  over  a range  of  temperatures. 

From  the  expression  for  the  width-and-shift  operator  in  Chapter  3,  it  is  clear  that 
this  operator  is  nearly  linear  in  the  electron  density  ne.  The  Debye  screening  cutoff 
introduces  a slight  deviation  from  this  linear  dependence.  Thus,  the  width  and  shift 
increases  with  increasing  electron  density.  The  physical  reasoning  for  this  behavior 
is  that  as  the  electron  density  increases,  a larger  number  of  binary  collisions  occur 
between  the  radiator  and  perturbing  electrons  during  the  lifetime  of  the  upper  state  of 
the  transition,  thereby  resulting  in  larger  electron-broadening  effects.  This  behavior  is 
shown  in  Figure  4.2  using  a full-Coulomb  treatment  of  the  perturbing  electron- radiator 
interaction  and  in  Figure  4.3  using  a dipole  approximation  to  this  interaction.  Therefore, 
we  have  shown,  as  was  stated  at  the  end  of  Section  4.1,  that  the  line  profiles  in  the  full- 
Coulomb  calculation  experience  a significant  amount  of  shift  and  distortion  while  the 
line  profiles  in  the  dipole-approximation  calculation  experience  very  little  shift  and  are 
relatively  symmetric.  When  the  density  becomes  sufficiently  high,  the  shift  and  possibly 
the  distortion  in  the  line  profile  should  be  noticeable  in  the  experimental  data,  and  thus 
a dipole-approximation  calculation  is  no  longer  adequate. 

The  large  shift  in  the  full-Coulomb  calculation  of  the  spectral  lines,  in  addition  to 
being  a useful  density  diagnostic,  gives  rise  to  some  very  interesting  physics.  Because 
a higher-lying  spectral  line  in  a Rydberg  series  shifts  by  a larger  amount  than  the 
adjacent  lower  spectral  line  of  the  series,  they  can  begin  to  merge  together  when  the 
density  becomes  large  enough.  Further  details  on  the  process  of  line  merging  are  given 


in  the  next  section. 
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Temeperature  (eV) 


Figure  4.1:  This  figure  shows  level  widths  and  shifts  due  to  electron  broadening  for  the 
Ar+17  Lyman  series  as  a function  of  temperature  at  an  electron  density  of  1 x 1024  cm-3. 
Note  that  the  values  of  the  widths  and  shifts  increase  with  principal  quantum  number  and 
decrease  with  increasing  orbital  angular  momentum  quantum  number,  a)  The  widths 
are  shown,  b)  The  shifts  are  shown. 
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Figure  4.2:  Line  shapes  of  hydrogen-like  argon  utilizing  a full-Coulomb  calculation  of 
the  electron-radiator  interaction  are  shown  as  a function  of  electron  density  with  a 
temperature  of  k qT  = 1 keV.  The  width  and  shift  of  the  line  shapes  increases  as  density 
increases.  Also,  note  that  the  line  shapes  contain  a distortion  due  to  the  differing  shifts 
of  the  levels,  a)  The  Ar+17  Ly -(3  line  is  shown,  b)  The  Ar+17  Ly-y  line  is  shown. 
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Figure  4.3:  Line  shapes  of  hydrogen-like  argon  calculated  in  the  dipole  approximation 
of  the  electron-radiator  interaction  are  shown  as  a function  of  electron  density  with  a 
temperature  of  k qT  = 1 keV.  Note  that  there  are  no  shifts  and  small  asymmetries  in  the 
line  shapes  under  the  dipole  approximation,  a)  The  Ar+17  Ly -/?  line  is  shown;  b)  The 
Ar+17  Ly-7  line  is  shown. 
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Although  strongly  dependent  on  the  electron  density,  it  can  be  seen  from  Figure  4.1 
that  the  width-and-shift  values  due  to  electron  broadening  are  relatively  weakly  depen- 
dent on  the  electron  temperature.  This  figure  shows  that  as  the  temperature  of  the 
electrons  increases,  the  width-and-shift  values  decrease.  The  physical  reasoning  for  this 
is  relatively  simple.  As  the  plasma  electrons  increase  in  temperature,  the  average  speed 
of  the  electrons  in  the  plasma  also  increases.  Therefore,  as  the  electrons  speed  past  the 
radiator  in  their  hyperbolic  trajectory,  there  is  less  time  for  them  to  interact  with  the 
radiator  thus  resulting  in  smaller  perturbation  effects  on  the  orbitals  of  the  radiator. 

What  is  the  effect  of  this  temperature  dependence  on  the  line  spectra?  The  answer  to 
this  question  depends  on  whether  a full-Coulomb  or  a dipole-approximation  treatment 
is  applied  to  the  perturbing  electron-radiator  interaction  in  the  calculation  of  the  width 
and  shift.  As  can  be  seen  in  Figure  4.4,  line  spectra  utilizing  a dipole  approximation 
change  by  small  amounts  with  changes  in  temperature.  However,  as  can  be  seen  in 
Figure  4.5,  line  asymmetries  in  the  theoretical  spectra  utilizing  a full-Coulomb  treatment 
of  the  electron-radiator  interaction  become  more  pronounced  as  the  plasma  electron 
temperature  decreases. 

These  asymmetries  arise  because  of  a significant  amount  of  differing  shifts  in  the 
angular  momentum  states  of  the  radiator  due  to  the  penetrating  monopole  contribution, 
as  can  be  observed  in  Figure  4.1.  The  theoretical  (3—  and  y-lines  in  Figure  4.5  represent 
3 p — Is  and  4 p — Is  transitions,  respectively.  Due  to  the  ion  Stark  splitting  of  the  line, 
the  upper  level  p states  are  mixed  with  other  angular  momentum  states,  thus  giving 
rise  to  differential  shifting  of  the  Stark  components.  This,  in  turn,  gives  rise  to  the 
asymmetries  seen  in  the  line  shape. 

4.3  Comparison  with  Other  Models 

In  this  section,  comparisons  are  made  between  the  results  of  the  SC  all-order 
model  and  a quantum  mechanical  model  to  second-order  in  the  radiator-perturbing 
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Figure  4.4:  Line  shapes  of  hydrogen-like  argon  calculated  in  the  dipole  approximation 
of  the  electron-radiator  interaction  are  shown  as  a function  of  temperature  with  a with 
a density  of  ne  = 1 x 1024  cm-3.  There  are  no  shifts  and  small  asymmetries  in  the  line 
shapes,  and  the  temperature  dependence  is  not  as  pronounced  as  is  seen  in  Figure  4.5. 
a)  The  Ar+17  Ly ~/3  line  is  shown;  b)  The  Ar+17  Ly-y  line  is  shown. 
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Figure  4.5:  Line  shapes  of  hydrogen-like  argon  utilizing  a full-Coulomb  calculation  of 
the  electron-radiator  interaction  are  shown  as  a function  of  temperature  with  a density 
of  ne  = 1 x 1024  cm-3.  The  differing  shifts  of  the  levels  cause  increasing  asymmetry  as 
the  temperature  decreases,  a)  The  Ar+17  Ly -/3  line  is  shown,  b)  The  Ar+17  Ly-y  line  is 
shown. 
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electron  interaction  (QM  second-order  model)  using  both  a full-Coulomb  and  a dipole- 
approximation  treatment  of  this  interaction.  Through  comparisons  of  the  results  of  the 
two  models,  we  demonstrate  the  plasma  conditions  under  which  the  second-order  model 
begins  to  break  down,  specifically  in  the  calculation  of  the  far  wings  of  the  line  shape 
and  calculations  of  adjacent  spectral  lines  in  a Rydberg  series  that  begin  to  overlap. 

We  begin  by  first  showing  the  comparisons  in  Figure  4.6  of  the  theoretical  Lyman-/? 
and  Lyman-y  line  profiles  from  both  QM  second-order  and  SC  all-order  models  using  a 
dipole-approximation  calculation  at  a typical  ICF  electron  density  and  temperature  of 
1 x 1024  cm-3  and  1 keV,  respectively.  We  can  see  that  the  QM  second-order  model  gives 
a line  profile  which  is  noticeably  more  broad  than  our  SC  all-order  model.  It  has  also 
been  shown  that  in  these  QM  second-order  models,  a dipole-approximation  calculation 
generally  overestimates  the  width  of  a line  [18].  As  can  be  seen  in  Figures  4.2  and  4.3, 
there  does  not  seem  to  be  quite  as  large  a discrepancy  in  the  SC  all-order  model.  In 
fact,  the  comparisons  of  line  shapes  between  versions  of  the  SC  all-order  model  using 
dipole-approximation  and  full-Coulomb  treatments  of  the  electron-radiator  interaction 
show  that  the  dipole-approximation  line  shapes  are  somewhat  more  narrow  than  the 
full-Coulomb  line  shapes.  In  Figure  4.7,  we  also  show  comparisons  of  full-Coulomb 
theoretical  Lyman-/?  and  Lyman-y  line  profiles  from  both  the  SC  all-order  and  the  QM 
second-order  models  at  the  same  density  and  temperature  stated  above.  We  see  that 
these  two  line  profiles  are  in  good  agreement  having  basically  the  same  width  and  shift. 

As  plasmas  become  more  strongly  coupled,  as  is  the  case  when  the  temperature 
decreases  or  the  density  increases,  the  validity  of  a second-order  theory  comes  into 
question.  If  we  keep  the  temperature  at  1 keV  and  raise  the  electron  density  to 
5 x 1024  cm-3,  the  comparison  of  the  full-Coulomb  SC  all-order  and  QM  second-order 
models  in  Figure  4.8  still  show  good  agreement  for  the  isolated  line  shape  calculation. 
However,  this  does  not  hold  true  when  mixing  is  included  between  the  upper-state  man- 
ifolds of  adjacent  lines,  as  is  demonstrated  at  the  end  of  this  section.  If  we  instead  keep 
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Figure  4.6:  This  figure  shows  comparisons  of  hydrogen-like  argon  spectra  from  the  dipole- 
approximation  version  of  the  SC  all-order  and  QM  second-order  models  at  an  electron 
density  and  temperature  of  1 x 1024  cm-3  and  1 keV,  respectively.  Note  that  in  the  dipole 
approximation,  the  QM  second-order  model  calculates  line  shapes  which  are  broader 
than  those  from  the  SC  all-order  model,  a)  This  plot  shows  the  Lyman-/3  line  of  Ar+17. 
b)  This  plot  shows  the  Lyman-7  line  of  Ar+17. 
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Figure  4.7:  This  figure  shows  comparisons  of  hydrogen-like  argon  spectra  from  the  full- 
Coulomb  SC  and  QM  models  at  an  electron  density  and  temperature  of  1 x 1024  cm-3 
and  1 keV,  respectively.  Note  that  the  two  models  are  in  good  agreement  for  these  line 
shapes,  a)  This  plot  shows  the  Lyman-/?  line  of  Ar+17.  b)  This  plot  shows  the  Lyman-7 
line  of  Ar+17. 
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the  density  at  1 x 1024  cm-3  and  lower  the  temperature  to  600  eV,  we  see  from  Figure  4.9 
that  the  full-Coulomb  SC  and  QM  models  are  still  in  good  agreement.  However,  their 
are  some  minor  differences  in  the  width  and  shift  of  the  lines.  In  fact,  we  do  not  see 
any  significant  differences  between  the  models  until  the  temperature  drops  to  or  below 
400  eV,  as  is  seen  in  Figure  4.10  and  the  corresponding  semilog  plot  in  Figure  4.11. 
Note  that  the  wings  of  the  line  are  beginning  to  differ  significantly  in  the  semilog  figure. 
According  to  the  conditions  set  forth  by  Cooper,  Kelleher,  and  Lee  [45] 

kBr  > Z-^~  ^Ryd  (4.1) 

nz 

the  temperature  limit  for  a second-order  calculation  is  460  eV  for  a Lyman-/?  line  and 
260  eV  for  a Lyman-y  line.  This  condition  states  that  the  momentum  of  the  perturbing 
electron  must  be  much  larger  than  any  change  in  its  momentum.  If  this  condition 
is  violated,  it  basically  means  that  strong  collisions  are  becoming  important  in  the 
determination  of  the  line  profile. 

For  perturbing  electrons  which  penetrate  the  orbitals  of  the  radiator  and  result  in 
large  shifts  and  distortions  in  the  line  profiles,  we  see  that  the  SC  all-order  model  and 
the  QM  second-order  model  still  tend  to  agree  with  one  another  in  the  amount  of  shift 
and  distortion  over  a wide  range  of  temperatures.  This  agreement  might  seem  odd 
because  electron  collisions  which  penetrate  the  orbital  may  seem  to  be  strong.  However, 
the  strength  of  an  electron  collision  depends  on  the  duration  of  the  collision  and  the 
electron’s  proximity  to  the  radiator.  The  duration  of  the  collision  is  an  important  factor 
because  the  electron  speeds  up  as  it  passes  by  the  radiator,  thus  decreasing  the  amount 
of  time  during  which  it  interacts  with  the  radiator.  Therefore,  fast  electrons  which 
penetrate  the  orbitals  can  still  experience  weak  collisions,  and  thus  the  center  region  of 
the  line  shapes  should  be  similar. 
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Figure  4.8:  This  figure  shows  comparisons  of  hydrogen-like  argon  spectra  from  the  full- 
Coulomb  SC  all-order  and  QM  second-order  models  at  an  electron  density  and  temper- 
ature of  5 x 1024  cm-3  and  1 keV,  respectively.  Note  that  the  two  models  are  still  in 
good  agreement  at  this  density,  a)  This  plot  shows  the  Lyman-/?  line  of  Ar+17.  b)  This 
plot  shows  the  Lyman-7  fine  °f  Ar+17. 
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Figure  4.9:  This  figure  shows  comparisons  of  hydrogen-like  argon  spectra  from  the  full- 
Coulomb  SC  all-order  and  QM  second-order  models  at  an  electron  density  and  temper- 
ature of  1 x 1024  cm-3  and  600  eV,  respectively.  Note  that  the  two  models  are  still  in 
good  agreement  at  this  temperature,  a)  This  plot  shows  the  Lyman-/?  line  of  Ar+17.  b) 
This  plot  shows  the  Lyman-7  line  of  Ar+17. 
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Figure  4.10:  This  figure  shows  comparisons  of  hydrogen- like  argon  spectra  from  the 
full-Coulomb  SC  all-order  and  QM  second-order  models  at  an  electron  density  and  tem- 
perature of  1 x 1024  cm-3  and  400  eV,  respectively.  Note  that  the  two  models  are 
beginning  to  diverge  in  accordance  with  the  temperature  limit,  a)  This  plot  shows  the 
Lyman-/?  line  of  Ar+17.  b)  This  plot  shows  the  Lyman-7  line  of  Ar+17. 
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Figure  4.11:  This  figure  shows  comparisons  of  hydrogen-like  argon  spectra  from  the  full- 
Coulomb  SC  all-order  and  QM  second-order  models  on  a semi-log  plot  at  an  electron 
density  and  temperature  of  1 x 1024  cm-3  and  400  eV,  respectively.  Note  that  the  two 
models  are  beginning  to  diverge  in  the  wings  when  A ui  > 2u Jpiasma-  Also,  the  energy 
scales  are  not  the  same  on  the  two  plots  for  reasons  of  comparison,  a)  This  plot  shows 
the  Lyman ~/3  line  of  Ar+17.  b)  This  plot  shows  the  Lyman-y  line  of  Ar+17. 
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We  have  shown  that  the  results  of  the  all-order  theory  and  the  second-order  theory 
are  generally  in  agreement.  Now,  we  consider  cases  in  which  the  SC  all-order  model  is 
needed  because  the  results  of  this  model  begin  to  show  deviations  from  the  results  of  a 
QM  second-order  model.  Remember  that  in  addition  to  the  temperature  limit  imposed 
on  the  second-order  model  in  Equation  4.1,  the  second-order  model  also  begins  to  break 
down  when  Auj  > 2u Jpiasma.  For  example,  look  at  the  semilog  plots  in  Figure  4.11. 
Because  ujpiasrna  ~ 37  eV  at  a density  of  1 x 1024  cm-3,  we  see  that  most  of  the  difference 
occurs  out  in  the  wings  of  the  line  shape  beyond  2u)p[asma  from  line  center.  Also,  low- 
Z radiators,  such  as  hydrogenic  helium,  can  emit  spectral  lines,  for  example  Balmer 
lines,  that  have  overall  widths  much  larger  than  twice  the  plasma  frequency.  This  limit 
can  also  be  violated  when  attempting  to  calculate  adjacent  lines  of  a Rydberg  series 
simultaneously.  This  type  of  calculation  is  required  when  considering  the  merging  of 
two  of  more  spectral  lines. 

As  the  width  and  especially  the  shift  of  the  spectra  in  a Rydberg  series  increases 
with  both  increasing  electron  density  and  increasing  principal  quantum  number,  the 
nth  member  of  this  series  shifts  toward  the  (n  — l)th  member.  When  the  electron  den- 
sity becomes  large  enough,  the  Stark  splitting  of  adjacent  energy  levels  due  to  the  ion 
microfield,  together  with  the  plasma  induced  shift  in  the  lines,  cause  these  series  mem- 
bers to  begin  merging.  As  this  happens,  the  theoretical  spectra  for  both  the  nth  and 
(n  — l)th  members  of  the  Rydberg  series  can  no  longer  be  calculated  independently  of 
one  another.  Off-diagonal  matrix  elements  between  the  upper  states  of  the  nth  and 
(n  — l)th  spectral  series  members  are  on  the  order  of  the  diagonal  matrix  elements,  and 
mixing  between  the  upper  states  of  these  series  members  is  significant.  For  example, 
the  magnitude  of  the  matrix  element  Bip^p  is  about  half  that  of  Bip^p  . Because  of  this 
mixing,  these  series  members  merge  together  to  form  a composite  spectral  feature  which 
does  not  obey  a linear-shift  model. 
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As  evidence  of  the  effect  of  mixing  between  the  upper  states  of  the  Lyman  /?-,  7-,  and 
Wines  of  Ar+17,  we  show  in  Figure  4.12  a comparison  of  this  composite  feature  with  and 
without  mixing  between  the  upper-state  manifolds  using  the  SC  all-order  model  with 
a full-Coulomb  electron-radiator  interaction  at  an  electron  density  and  temperature  of 
5 x 1024cm-3  and  1 keV,  respectively.  The  calculation  without  mixing  is  often  referred 
to  as  the  linear  model  because  we  effectively  are  adding  the  intensities  of  the  individual 
series  members  in  which  the  density  dependence  of  the  shift  of  these  members  are  nearly 
linear.  A drawback  of  the  linear  model  is  that  the  shift  of  a Rydberg  series  member  to 
lower  energy  increases  with  increasing  principal  quantum  number  and  increases  nearly 
linear  with  density,  thus  providing  circumstances  in  which,  for  example,  a (5-line  can 
shift  past  a 7-line.  The  simultaneous  calculation  of  the  the  Lyman  series  members  with 
the  inclusion  of  mixing  destroys  this  linear  shift  dependence  on  the  electron  density,  and 
thus  is  often  referred  to  as  the  nonlinear  model.  The  mixing  between  the  upper-state 
manifolds  of  the  Rydberg  spectral  series  members  in  the  nonlinear  model  prevents  one 
member  of  a spectral  series  from  shifting  past  another  member. 

Finally,  let  us  consider  differences  between  the  SC  all-order  model  and  the  QM 
second-order  model  using  a full-Coulomb  treatment  of  the  electron-radiator  interaction. 
Figures  4.13  and  4.14  show  comparisons  of  the  composite  Lyman-/?7(5  feature  at  the 
three  different  electron  densities  of  1 x 1024  cm-3  , 2.5  x 1024  cm-3,  and  5 x 1024  cm-3 
with  an  electron  temperature  of  1 keV.  As  a reference,  we  also  include  the  calculation 
of  the  linear  SC  all-order  model  in  the  figures.  We  see  that  as  the  electron  density 
increases,  the  observable  difference  on  the  blue  side  of  this  composite  feature  increases 
and  is  large  enough  to  be  experimentally  noticeable  at  densities  near  5 x 1024  cm-3.  The 
difference  shows  that  the  SC  all-order  model  predicts  a process  of  merging  which  occurs 
more  slowly  with  increasing  electron  density.  The  reason  for  this  lies  in  the  frequency- 
dependent  resolvent  in  Equation  B.8  of  Appendix  B.  In  the  QM  second-order  theory 
mentioned  in  this  work,  the  electron-radiator  interaction  term  present  in  the  denominator 
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feature  calculated  with  the  full-Coulomb  SC  all-order  model  with  and  without  mix- 
ing between  the  upper-state  manifolds  of  the  individual  lines  at  an  electron  density  of 
5 x 1024  cm-3  and  an  electron  temperature  of  1 keV. 

is  replaced  by  the  first-order  static  shift.  The  reason  for  using  this  replacement  instead  of 
dropping  the  interaction  term  is  that  the  value  of  this  term  approximates  the  contribution 
from  the  monopole  part  of  V^r.  This  term  can  be  of  a sufficiently  large  value  such 
that  it  cannot  be  treated  as  a perturbation.  Therefore,  the  QM  second-order  theory 
depends  on  a resolvent,  and  thus  a time  development  operator,  which  is  independent 
of  the  dynamic  nature  of  the  perturbing  electron-radiator  interaction  V^\.  This,  in 
turn,  underestimates  the  dynamic  behavior  of  the  mixing  between  the  upper  states  and 
overestimates  the  rate  at  which  the  members  of  the  series  merge  together  as  the  electron 
density  increases. 
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Figure  4.13:  This  figure  shows  the  comparison  of  the  theoretical  composite  Lyman- 0^6 
feature  of  Ar+17  using  the  nonlinear,  full-Coulomb  SC  all-order  model  (includes  mixing), 
the  nonlinear,  full-Coulomb  QM  second-order  model  (includes  mixing),  and  the  full- 
Coulomb  linear  calculation  (no  mixing)  which  is  nearly  identical  in  both  QM  and  SC 
models  at  a temperature  of  1 keV  over  the  density  range  shown,  a)  The  spectra  are 
calculated  at  an  electron  density  of  1 x 1024  cm-3,  b)  The  spectra  are  calculated  at  an 
electron  density  of  2.5  x 1024  cm-3. 
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Figure  4.14:  This  figure  shows  the  comparison  of  the  theoretical  composite  Lyman- (3j5 
feature  of  Ar+17  using  the  nonlinear,  full-Coulomb  SC  all-order  model  (includes  mixing), 
the  nonlinear,  full-Coulomb  QM  second-order  model  (includes  mixing),  and  the  full- 
Coulomb  linear  calculation  (no  mixing)  which  is  nearly  identical  in  both  QM  and  SC 
models  at  an  electron  density  and  temperature  of  1 x 1024  cm-3  and  1 keV,  respectively. 


CHAPTER  5 
CONCLUSIONS 


In  this  dissertation,  we  presented  a semiclassical,  all-order,  full-Coulomb  electron- 
broadening model  which  includes  the  effects  of  strong  collisions  through  the  inclusion  of 
the  terms  to  higher-order  in  the  electron-radiator  interaction.  To  calculate  theoretical 
spectra  for  higher-Z  radiators  and  to  include  contributions  from  strong  collisions,  we 
developed  a model  that  differs  from  the  VCS  charged  radiator  model  in  these  important 
ways: 

• The  interaction  between  the  perturbing  electron  and  the  radiator  is  calculated 
using  a multipole,  or  full-Coulomb  calculation. 

• Relativistic  wavefunctions  and  the  corresponding  energy  level  structure  from  a 
relativistic  atomic  physics  package  are  employed  to  describe  the  atomic  physics  of 
higher-Z  radiators. 

• A minimum  separation  cutoff  is  incorporated  into  the  model.  This  cutoff  is  the 
Landau  length  modified  by  the  radiator  charge  and  quantum  diffraction  effects. 
The  cutoff  corrects  for  the  singularity  caused  by  the  treatment  of  the  electrons  as 
classical  particles. 

• The  distribution  function  is  generalized  to  provide  temperature-dependent  changes 
in  the  relative  populations  of  the  radiator  states  modified  by  the  interaction 
between  the  perturbing  electron  and  the  radiator. 

• The  no-quenching  approximation,  which  prohibits  nonradiative  transitions  caused 
by  perturbations  between  states  of  differing  principle  quantum  number  n,  must  be 
relaxed  to  allow  for  nonradiative  transitions  between  upper  states  of  spectral  lines 
which  are  beginning  to  merge  together. 

The  results  of  this  model  were  compared  with  a quantum  mechanical,  second-order,  full- 
Coulomb  electron-broadening  model  which  does  not  treat  strong  collisions  as  thoroughly 
as  our  model  does  due  to  the  exclusion  of  the  higher-order  terms.  We  considered  the 
impact  of  strong-collision  effects  through  the  higher-order  terms  on  line  merging  and 
other  related  physical  phenomenon. 

We  discovered  that  the  inclusion  of  strong-collision  effects  through  the  higher-order 
terms  in  the  interaction  can  have  a significant  effect  on  the  frequency  dependence  of 
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the  width  calculation  in  the  wings  of  a line  profile  emitted  by  highly  ionized,  mid-Z 
radiators.  It  also  has  a significant  effect  on  how  quickly  line  profiles  merge  together  as 
the  density  increases.  As  was  shown,  the  spectral  lines  from  the  all-order  model  tend  to 
merge  more  slowly  than  the  line  profiles  from  the  second-order  model. 

As  electron  densities  in  inertial  confinement  fusion  plasmas  increase,  this  model  has 
the  potential  to  be  useful  as  a tool  in  plasma  spectroscopy  diagnostics.  As  densities 
become  sufficiently  high,  line  profiles  begin  to  merge  together,  and  the  shape  of  the  com- 
posite feature  seems  to  have  a significant  dependence  on  strong-collision  effects  through 
the  inclusion  of  the  higher-order  terms.  Also,  the  dipole  approximation  has  been  shown 
to  be  inadequate  at  higher  electron  densities  where  line  shifts  and  line  merging  are  impor- 
tant. This  makes  sense,  because  most  of  the  strong-collision  effects  occur  through  the 
monopole  interaction  of  the  electron-radiator  interaction  when  the  perturbing  electron 
penetrates  the  orbitals  of  a mid-Z  radiator.  These  results  are  significant  because  plasma 
diagnostics  may  be  extended  to  higher  densities  as  researchers  work  to  achieve  “break 
even”  energy  creation  in  fusion  research. 

Differences  in  merged  spectral  lines  due  to  strong  collisions  have  the  possibility 
of  altering  what  we  perceive  as  the  distinction  between  continuum  and  bound  states. 
Because  ionization  balance  is  sensitive  to  any  changes  in  the  energy  level  structure  [28-30] 
and  opacity  is  sensitive  to  the  energy  location  and  intensity  of  line  emission,  an  under- 
standing of  the  effects  of  strong  collisions  on  the  merged  composite  spectral  features 
which  lie  near  the  continuum  edge  can  be  important. 

Even  at  lower  densities,  we  have  shown  that  the  frequency  dependence  of  the  width 
in  the  wings  of  a line  profile  are  significantly  dependent  on  the  effects  of  strong  collisions 
through  the  higher-order  terms.  This  can  have  significant  effects  on  Rosseland  group 
opacities  [31].  They  are  very  sensitive  to  changes  in  the  intensity  of  the  line  profile  in 
the  wings  because  they  are  proportional  to  a harmonic  mean  of  the  absorption  coeffi- 
cient. Rosseland  group  opacities,  in  turn,  are  important  in  the  calculation  of  multi-group 
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diffusion  radiative  transfer  effects.  These  effects,  which  are  important  in  determining  the 
amount  of  radiation  that  passes  through  different  portions  of  the  frequency  spectrum, 
can  be  of  importance  in  both  astrophysical  studies  and  in  fusion  research. 

There  are  areas  that  need  further  examination  in  this  model.  First,  as  the  density  in 
the  plasma  increases,  the  one-electron  approximation,  where  we  allow  only  one  strong 
electron  collision  with  the  radiator  at  a time  along  with  the  treatment  of  electron  cor- 
relations through  a static  Debye  shielding  model,  become  questionable  and  therefore 
requires  a more  rigorous  treatment.  Second,  in  order  to  calculate  spectral  lines  which 
have  more  than  one  final  state,  the  model  must  be  generalized  to  include  lower-state 
broadening  in  addition  to  upper-state  broadening.  Third,  the  plasma  electrons  do  not 
interact  with  the  radiator  in  isolation  with  the  presence  of  nearby  ions,  thus  requiring 
at  least  a calculation  of  the  width-and-shift  operator  in  the  presence  of  the  ion  electric 
microfield.  Fourth,  an  examination  of  the  quantum  effects  of  intermediate  momentum 
states  of  the  perturbing  electrons  on  line  merging  should  be  performed.  Finally,  at  this 
time,  we  do  not  have  conclusive  experimental  evidence  to  confirm  the  behavior  of  line 
merging  as  predicted  by  the  theoretical  models  as  of  yet.  Our  NLUF  (National  Laser 
User  Facility)  campaign  is  currently  underway  on  the  OMEGA  laser  with  the  intention  of 
generating  sufficiently  high  densities  and  temperatures  by  varying  laser  pulse  shapes,  fill 
pressures  of  the  microballoon,  and  the  thickness  of  the  ablator.  Preliminary  results  indi- 
cate electron  densities  somewhat  above  2 x 1024  cm-3  and  temperatures  above  1.5keV, 
and  some  very  preliminary  fits  suggest  that  the  nonlinear  models  fit  somewhat  better 
than  the  linear  model.  Also,  with  the  startup  of  the  NIF  (National  Ignition  Facility) 
laser  in  2003,  the  increased  laser  energy  on  target  should  allow  us  to  probe  plasmas 
in  higher  temperature  and  density  regimes  and  allow  us  to  more  conclusively  test  the 
models  in  question. 


APPENDIX  A 

EQUATION  OF  MOTION  FOR  D{u,  e) 


In  this  appendix,  we  use  Zwanzig’s  projection  operator  formalism  [46]  to  derive  the 
equation  of  motion  for  D(lo,  e).  We  then  interpret  the  physical  meaning  of  each  of  the 
terms  in  the  equation  of  motion.  From  these  terms,  we  then  develop  a model  for  both 
ion  broadening  and  electron  broadening  of  the  spectral  line  shape. 

To  calculate  the  Stark-broadened  line  shape,  we  need  to  derive  the  equation  of  motion 
for  Equation  2.22 

D(u,  e)  = — J d t S"-*#  Trp  I <5  (e  - Esi)  p'e~iH  ^hdeiH  */*]  . (A.l) 

To  simplify  the  form  of  this  expression  and  the  corresponding  equation  of  motion,  we 
apply  the  properties  of  a Liouville  operator 


Lf  = ~n 


H ./ 


t/h  f gi H t/h 


(A.2) 


and  perform  the  Laplace  transform.  Thus,  D(u,e)  becomes 

<A,3) 

In  the  application  of  the  Zwanzig  projection  operator  method,  it  useful  to  define 

p (e)  ^ 5(€-e")p  (a-4) 

that  is  interpreted  as  the  equilibrium  state  distribution  function  for  a radiator  immersed 
in  a plasma  constrained  to  have  a microfield  of  the  value  e.  It  is  also  useful  to  define 
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the  Laplace  transform  of  the  time  evolution  of  the  radiator  dipole  as 

X(u)=T±—d.  (A. 5) 

(W  - L) 

so  that  D(u,e)  now  has  the  form 

£>(^,e)  = Trp{p'(e)X(a;)}  . (A.6) 

Utilizing  the  Zwanzig  projection  operator  method,  we  can  define  the  projection  oper- 
ator for  the  line  broadening  problem  as  having  the  action 

PY  = J de  6 (e-  Esi)  p~\  Trp  {/(e)  Y } (A.7) 

where 

Pe,r  = Trpp'(e)  (A. 8) 

is  the  radiator  reduced  distribution  function  in  the  presence  of  the  ion  microfield.  This 
projection  operator  selects  the  portion  of  the  operator  Y in  the  radiator-plasma  system 
contributing  to  the  ensemble  average  over  the  possible  values  of  the  ion  microfield. 
Applying  this  projection  operator  to  the  operator  X(oS)  as  defined  above,  we  have 

PX(uj)  = J de  5 (e  - Esi)  P~\  D(",  <0  • (A.9) 

Because  P is  a projection  operator,  it  satisfies  the  relation  P 2 = P.  From  the  form 
of  P,  it  is  also  easy  to  show  that  Pd  = d because  d does  not  depend  on  perturber 
coordinates.  A projection  operator  of  the  form  Q = 1 — P picks  out  the  portion  of  a 
system  operator  which  is  excluded  when  applying  the  projection  operator  P.  We  now  use 
these  two  projection  operators  to  derive  the  equation  of  motion  for  D{ui,e).  Applying 
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P to  Equation  A. 5 and  then  Q to  Equation  A. 5,  we  get  a pair  of  equations 


(u>  - PLP)  PX  - ( PLQ ) QX  = d (A.  10) 

(w  - QLQ)  QX  - (QLP)  PX  = 0 (A.ll) 

By  solving  these  two  equations  simultaneously,  we  find  the  equation  of  motion  for  PX 

(c j - PLP  - PLQ  {u  - QLQy1  QLP ) PX  = d.  (A.  12) 


Finally,  by  inserting  the  definition  of  the  projection  operator  P into  the  equation  of 
motion  for  PX , we  arrive  at  the  formal  equation  of  motion  for  D(uj,e) 

uD(u,  e)  — J de  |B(e,  e)  + M(e,  e;n;)|  D(u,  e)  = pe,rd  (A. 13) 

where  the  instantaneous  or  static  contributions  to  the  broadening  of  the  line  shape  are 
contained  in  the  operator 

B(e,e')  = Trp  |p'(e)L<5(e'  - Esi)^  p~Jr  (A.14) 

and  the  dynamic  perturber  collision  contributions  to  the  broadening  of  the  line  shape 
are  contained  in  the  operator 

M(e,  e-u)  = Trp  jp'(e)LQ-_-L_  QL6(e'  - £si)}  p~}r  . (A.15) 

The  dynamic  broadening  operator  M(e,e  ;tj)  is  very  complicated.  Therefore, 
to  illustrate  the  method  at  hand,  we  begin  by  analyzing  the  contributions  to  the 
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instantaneous  broadening  operator  B(e,e').  Inserting  the  definition  of  p'(e),  we  have 


')  = ^yTrp| 6{e  - Eai)p' L5(e'  - Esi)}  pj7 


This  expression  can  be  rewritten  in  the  form 


(A.  16) 


£(«.«')  = ;^Trp{L<S(e-.E,i)p'ii(e'-.E,l)}p<,1r 

~W)Tl”  {|W(e  “ E’,)]p'S(e' " p~lr  (A17> 

because  L does  not  commute  with  S(e  — Esi).  Substituting  in  the  terms  of  the  Liouville 
operator,  this  expression  can  be  shown  to  have  the  form 


B(e,  e)  = 6(e  — e) 


Lr  + W)TrP  {Lj5(e  - Esi )}  pj 


(A.18) 


The  contribution  to  L;  from  the  dipole  interaction  for  both  electrons  and  ions  can  be 
separated  off  and  together  give  the  Stark  broadening  operator,  where  Esi  is  identified  as 
the  ion  field  screened  by  the  electrons  [47].  Therefore, 


B(e,  e ) — S (e  — e')  [ Lr  + Lj)r  (e)  + Be(e)] 


(A.19) 


where 

B.w  = ^Trp  _ fy/}  p-1  (A. 20) 

is  the  part  of  the  instantaneous  broadening  operator  containing  the  mean  field  effects 
of  the  perturbing  electrons  on  the  radiator.  The  primes  on  LKe,r  correspond  to  the 
electron-radiator  interaction  minus  the  monopole  and  dipole  contributions.  Thus,  the 
equation  of  motion,  which  is  still  exact,  now  has  the  form 


[u>  — Lr  — Litr(e)  — Be(e)]D(uj,e)  — / de,M(e,  e \ lu)D(u,  e')  = p£,rd . (A. 21) 
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As  was  noted  above,  the  analysis  of  is  more  complex  and  can  be  done 


qualitatively  as  follows.  First,  note  that  the  appearance  of  QL  assures  that  only  the 
interaction  parts  of  L occur  (i.e.,  QL  — ► QL/).  Thus,  M(e,e';u;)  has  contributions 
describing  broadening  by  ion  dynamics,  electron  dynamics,  and  cross  terms  representing 
their  correlation.  As  a first  approximation,  we  assume  that  there  are  no  dynamical 
correlations  between  the  electrons  and  the  ions.  Therefore,  it  can  be  shown  that 


- J de' Mt(e,  e';u)D(u,  e)  = p^rd  (A. 23) 

which  is  quoted  and  discussed  in  the  text,  noting  that  the  various  broadening  mechanisms 
have  been  isolated  and  identified. 

Using  an  ion  motion  assumption  in  which  the  transitions  among  different  field  values 
are  instantaneous  and  uncorrelated  (Poisson  process),  the  ion-dynamics  collision  term 
can  be  written  in  the  form  [38] 


M(e,  e ; u)  — ■>  Mj(c , e']u>)  + 6(e  — e')Me(uj , e) 


(A. 22) 


where  M*(e,  and  Me(e;  u>)  describe  the  dynamical  effects  of  the  ions  and  electrons, 
respectively.  Using  this  approximation,  the  equation  of  motion  has  the  form 


[a )-Lr-  Li>r(e)  - Be(e)  - Me(uj , e)]  D(u>,  e) 


(A. 24) 


where  u is  directly  proportional  to  the  momentum  autocorrelation  function  of  the  ions 
and  is  a measure  of  the  importance  of  ion  motion  effects  on  the  line  shape.  With  this 
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definition,  we  see  that  pe<r  is  a stationary  point  of  Mt.  or 


/ 


d e'  Mi{e,e'-uj)p£'  = 0, 


(A.25) 


Now,  if  we  define 


G(uj,  e)  — 


• 1 

uj  — Lr  — Litr(e)  — Be(e)  — Me(u,  e)  — iv  ’ 


(A. 26) 


the  equation  of  motion  for  -D(w,  e)  has  the  form 


D(uj,  e)  = G(u,  e) 


rd  — 'wpeirJ  d e' Q(e')pJ rD(uj,e') 


(A. 27) 


Next,  we  solve  the  equation  of  motion  and  find  an  expression  for  D(u,  e).  By  defining 
the  quantity  Y as 

Y = j de'  Q(e')p-}rD{u,e)  (A.28) 

and  substituting  Equation  A. 27  into  Y,  we  have 


Y = J d e'Q(e')pJrG(u,  e)  [/v  ,rd  - i "Pe',rY] 


(A. 29) 


Solving  this  expression  for  Y , we  find  that 


^ =l  + i,Jde-Q(e^^G(a,,eO/,.v/d6'Q(6'AvG("’e>-'^-  <A-30> 


By  noting  the  definition  of  Y and  substituting  this  expression  into  Equation  A. 27,  we 
find  that  after  some  simple  reduction 


D(u,e)  = G(uj1e)pe>7 


1 + ivf  de'Q(e')pJrG(uj,e')pe't 


d. 


(A.31) 
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Also,  from  the  definition  of  pe>  r,  it  can  be  shown  that 


G(u,e')  = pJrG{uj,e')pi 


€ ,r 


that  is  then  used  to  arrive  at  the  final  expression  for  D(u>,  e)  used  in  the  text, 
expression  has  the  form 


D(u,e)  = G(ijj,e)pe iT 


1 + iv  fd e' Q(e')G(uj,e') 


d. 


(A. 32) 
This  final 


(A. 33) 


APPENDIX  B 

ELECTRON-BROADENING  OPERATOR 


From  Appendix  A,  we  found  that  the  instantaneous  or  static  contribution  to  the 
electron-broadening  operator  has  the  form 


Bc(e)  = [j(e  - E«)p'] } 


P.'r 


(B.l) 


In  Appendix  A,  we  also  applied  the  approximation  that  there  are  no  dynamical  cor- 
relations between  the  ions  and  the  electrons.  Under  this  approximation,  the  dynamic 
portion  of  the  electron-broadening  operator  has  the  form 


Me(uj,  e)  = Tre 


1 

QL(e)Q 


(B.2) 


where 

Z/(e)  = Lr  + Le  + + L (B.3) 

describes  the  dynamics  of  the  electron-radiator  interactions.  The  presence  of  the  projec- 
tion operator  Q under  the  current  approximation  selects  the  part  of  the  occurrences  of  L 
in  the  numerator  containing  the  the  electron-radiator  interaction,  L^r  ■ Also,  Lr  is  the 
Liouville  operator  of  the  isolated  radiator,  Le  is  the  Liouville  operator  for  the  plasma 
electrons,  and  L,,r(e)  is  the  Liouville  operator  for  the  ion  Stark  splitting.  Note  that  the 
ion  degrees  of  freedom  have  been  integrated  out  because  they  have  been  assumed  to  be 
“static”  in  the  electron-broadening  operator. 

By  applying  a weak-coupling  approximation  between  the  ions  and  the  electrons  which 
is  given  by 

P = Pe,e,rPi  (B.4) 
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the  static  and  dynamic  parts  of  the  electron-broadening  operator  now  have  the  form 


Note  that  in  this  simplification,  we  have  used  the  definition  of  the  ion  microfield  proba- 
bility distribution  function,  and  peer  is  the  reduced  electron-radiator  density  matrix  in 
the  presence  of  the  ion  field  e. 

We  next  apply  a one-electron  approximation  in  which  we  neglect  electron  correlations 
and  simultaneous  strong  electron-radiator  interactions  by  two  or  more  electrons.  Under 
this  approximation,  the  occurrences  of  the  projection  operator  Q vanish,  and  the  static 
and  dynamic  parts  of  the  electron-broadening  operator  now  have  the  form 


To  compensate  for  the  removal  of  electron  correlations,  we  apply  a Debye  static  shielding 
approximation  to  the  calculation  of  the  electron-broadening  operator  [39]. 

Next,  we  reformulate  the  expression  for  the  dynamic  electron  collision  operator 
in  terms  of  a time  development  operator  using  the  interaction  representation.  In  the 
dynamic  electron-broadening  operator,  the  frequency  dependent  denominator  contains 
all  of  the  dynamic  or  time  dependent  properties  of  the  electron  collisions  with  the 


(B.5) 


and 


(B.6) 


B(e)  = n, Tii, 


(B.7) 


and 


M(lo,  e)  = ne Trle 


where 


L,U)  = L,+L^(t). 


(B.9) 
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radiator.  Through  an  inverse  Laplace  transform,  this  resolvent  term  can  be  expressed 
as 

- = -[[  dte^-Lr(e)-Lu-L[\\r+ie)t_  (B.10) 

UJ  — Lr(e)  — Lie  ~ Lie  r Jo 

However,  because  the  electron  broadening  of  the  line  shape  is  due  to  the  perturbing 
electron- radiator  interaction,  it  is  advantageous  to  rewrite  this  quantity  in  the  interaction 
representation.  Therefore,  the  resolvent  can  be  expressed  in  terms  of  a time  development 
operator 

1 r 00 

_ = -i/  dtei{u,-LrU-L'*+k)tU(t,0)  (B.ll) 

w-Lr(€)-Llt-L<»r  Jo 

where 

0)  = tfexp  |-i  J dfL^r(t)|  (B.12) 


and 


i(Lr (c)+Lie)i  t (1)  —i(Lr(e)+Lie)t 
^1  e,r  e 


(B  13) 


The  time  development  operator  and  its  dynamical  properties  are  defined  by  the  differ- 
ential equation 

i^T^  = £gr(*Mi,0),  (B.14) 

and  i9  represents  a time-ordered  treatment  of  the  exponential  similar  to  that  seen  in  a 
Dyson  expansion.  With  the  substitution  of  the  expression  for  the  resolvent,  the  dynamic 
electron-broadening  operator  can  be  expressed  as 


M(uj,  e) 


-i  ne 


l 


d£ ei(UJ-ir(€)+i£)t  Trle  ^L^r(t) 
X lA(t , 0)  ple,e,r  -^le^r  } Pe,r 


(B.15) 


where  the  invariance  of  the  trace  under  a time  translation  of  the  perturbing  electron 
coordinates  has  been  used  to  simplify  this  operator. 
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Because  we  are  interested  specifically  in  radiative  transitions  between  initial  and 
final  states,  we  now  apply  the  no-quenching  approximation.  From  Equation  2.27,  it 
is  observed  that  the  width-and-shift  operator  acts  on  the  radiator  dipole  moment  d. 
Under  the  no-quenching  approximation,  the  matrix  elements  of  the  perturbing  electron- 
radiator  interaction  between  initial  and  final  states  of  the  transition  are  zero.  This 
results  in  matrix  elements  only  between  initial  and  final  states  of  the  radiator  transition 
for  this  occurrence  of  the  radiator  dipole  moment,  noted  by  dj/.  If  we  operate  B(e) 
and  M(u),e)  on  d*/  and  expand  the  commutator  definitions  of  the  Liouville  operator, 
we  then  have 


5(e)djy  — ^eTf le-fUjg  r Pie,e,r  Pe 

7leTrie{pie]£  r pc  rdjy  Vje  r} 


(B.16) 


and 


(i) 

1 e.r 


/oo  ( 

At Trle  { Vffr(t)  U(t, 0) 

xpie,.,rVH\p-)d„ 

+ in,  / °°  Trie  {{/((,  0)  pu,.,r  V, 

xp.:'d 

r oo 

+ me  J d t Trie  | U(t,  0) 

X Pie,e,rPe,l  dif  V\l\  0)  } e'H*'rt 

r oo 

/ d Trle  {U(t,  0)  ple,e,r  p~x 
Jo 

xdlfv}»U<(t,0)V^T(t)} 


e'H€,rt 


(B.17) 


The  first  term  describes  the  broadening  on  the  upper  states  of  the  transition,  the 
last  term  describes  the  broadening  on  the  lower  states  of  the  transition,  and  the  two 
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intermediate  terms  in  M(lj,  e)  are  interference  terms  between  upper  and  lower  states. 
Because  of  the  no-quenching  approximation,  operators  to  the  left  of  d,/  act  on  the  initial- 
state  manifold  of  the  transition  and  operators  to  the  right  of  d*/  act  on  the  final-state 
manifold  of  the  transition. 

Because  this  dissertation  concentrates  on  line  spectra  which  have  no  significant  lower 
state  broadening,  we  apply  a no  lower-state  broadening  approximation  which  is  discussed 
in  Chapter  3.  Therefore,  we  concentrate  on  the  first  term  of  instantaneous  part  of  the 
electron-broadening  operator 


Btl"‘(e)d  = n.TMWe1,’.  p~'rd 


(B.18) 


and  the  first  term  of  the  dynamic  part  of  the  electron-broadening  operator 


r o° 

Mfirst(uj,  e)d  = — i ne  / d t «■+“)*  Trle  j V$r(t)  U(t , 0) 

X Ple,e,r  V$r  p£  d U^t , 0)  } jH , (B.19) 


where  the  time  development  operator  has  the  form 


U(t,0)  = i?exp  {-if  d t Vg}r(t') 


(B.20) 


and 


Vffr(t)  = e^+H^  V$r  e-W-'+HlJ 


(B.21) 


The  time  development  operator  satisfies  the  differential  equation 


ihEM.  = vmrm(tt0). 


(B.22) 


These  are  the  expressions  which  are  used  in  the  main  body  of  the  text  for  the  calculation 
of  the  line  shapes  of  interest  in  this  dissertation. 


APPENDIX  C 

DEUTSCH  POTENTIAL  AND  THE  QUANTUM  CUTOFF 


In  the  model  presented  in  this  dissertation,  we  consider  the  case  of  electrons  treated 
as  classical  particles  undergoing  collisions  with  a radiating  ion  which  gives  rise  to  broad- 
ening effects  on  the  emitted  spectral  lines.  However,  under  a classical  treatment  of 
the  radiator-electron  Coulomb  interactions,  the  potential  goes  to  negative  infinity  as 
the  distance  between  the  radiator  and  the  electrons  go  to  zero.  In  order  to  remove 
this  singularity,  a minimum  distance  cutoff  is  derived  by  applying  quantum  diffraction 
effects  through  the  calculation  of  the  pair  potential  introduced  by  Deutsch  [40].  In  this 
appendix,  we  sketch  out  the  form  of  this  pair  potential.  The  pair  potential  is  then 
used  to  define  a minimum  distance  cutoff  that  used  in  the  electron-broadening  model 
presented  in  this  dissertation. 

The  effective  pair  potential  Vy(r ) is  defined  as  a function  of  the  two-particle  radial 
distribution  function,  gij(r),  by 


For  the  radiator-electron  interactions,  the  radial  distribution  function  can  be  considered 
as  the  sum  of  contributions  from  bound  and  scattered  states  of  the  radiator 


Vij(r)  = — kBT  In  g^r) . 


(C.l) 


9ij{r)  = gs(r)  + gb(r). 


(C.2) 
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The  scattered  state  contribution  can  be  written  in  terms  of  Coulomb  wave  functions 
Fi(a,kr ) 


J0  1 = 0 


(21  + l)2 
/ + 1/2 


|Fi(q,  *t)|2 


(C.3) 


where  Ay  is  the  thermal  De  Broglie  wavelength 


\ ( * 
{MMT 


1/2 


(C.4) 


and 

—Ze2Ha 

a=  h2k 


(C.5) 


Z is  the  overall  charge  of  the  radiating  ion  and  /iy  is  the  reduced  mass  of  the  electron- 
radiator  pair,  which  is  effectively  equal  to  the  electron  mass  for  most  cases  of  interest. 
The  bound  state  contributions  can  be  written  in  terms  of  the  bound  state  wave  functions 
of  the  radiator  ^n;m(r , 9,  <j)) 


Qb(r)  = (27rA?  )3/2  JJ  sin  9 &9  dcf)  ^ \^nim{r,9,(f>)\2 . (C.6) 

n,  Z,  771 


We  now  define  the  minimum  distance  cutoff  used  in  the  calculation  of  the  electron- 
broadening operator  as 


(C.7) 


where  Vc(r)  is  the  bare  Coulomb  potential  between  the  radiator  and  the  electron.  Sub- 
stituting in  the  definitions  for  the  potentials,  we  have 


T min 


Z T ind 

In  gij(r  -*  0) 


(C.8) 


80 


where  r;ncj  is  the  Landau  length 

rind  = koT  ' ^C'9) 

The  use  of  Zrin<i  as  the  cutoff  would  effectively  limit  the  strength  of  the  perturbing 
electron-radiator  interaction  to  approximately 


V(1) 

1 e.r 


k bT 


(C.10) 


and  assures  that  the  strength  of  the  interaction  is  balanced  by  the  speed  of  the  electron 
as  it  passes  by  the  charged  radiator.  Also,  note  that  In  gij(r  — ► 0)  modifies  the  cutoff 
for  quantum  diffraction  effects  between  the  perturbing  electron  and  the  radiator  [40] 
because  there  is  overlap  of  the  wave  packet  of  the  perturbing  electron  with  the  radiator 
orbitals. 


APPENDIX  D 

IRREDUCIBLE  TENSOR  OPERATORS  AND  3N-J  SYMBOLS 


This  appendix  provides  a brief  overview  of  a formalism  that  is  very  useful  for  the 
treatment  of  angular  momentum  dependent  matrix  elements.  These  definitions,  iden- 
tities, and  sum  rules  are  be  employed  in  the  calculation  of  the  matrix  elements  of  the 
perturbing  electron-radiator  interaction.  First,  Wigner  3N-j  symbols  [43,48,49]  are  intro- 
duced in  terms  of  an  integral  over  spherical  harmonics, 


Y^YkqYi'm'SinOdQdcf)  = (-l)m 


k V 

q m' 


\ 
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There  are  several  useful  relations  and  symmetry  properties  that  are  employed  in  this 
work,  such  as 


/ . . 

3 2 Jl 


l m2 


m i 


. \ 

3 3 


m3 


( 


_1  V'l+J2+J3 


(-1) 


Jl 

rrii 


32 

m2 


. \ 
33 

m3  J 


(D.3) 


and 


/ . . . \ ( . ■ \ 


3\ 

32 

3 3 

_ ^_±yi+h+h 

3 1 

32 

3 3 

^ ~m i 

-m2 

~m3  j 

^ mi 

m2 

m3  y 

(D.4) 


81 


82 


The  following  relation  can  be  obtained  using  the  orthogonality  of  spherical  harmonic 
functions, 
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Also,  the  equation  rri\  + m2  + m3  = 0 must  be  satisfied  for  a 3-j  symbol  to  be  nonzero. 
He  are  several  useful  sum  rules: 
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where  S(j\j2j3)  is  +1  if  the  triangle  relations  are  satisfied. 

The  Wigner-Eckart  Theorem  [41,42]  is  employed  to  express  the  matrix  elements  of 
our  interaction  potential  in  terms  of  reduced  matrix  elements  as  follows: 
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where  Tq ^ is  an  irreducible  tensor  operator,  j and  / are  the  angular  momentum  indices 
of  the  coordinates  relevant  to  the  operator,  and  a and  a'  represent  all  the  other  quantities 
(energy,  spin,  etc.)  that  define  their  respective  states. 

A particular  example  of  reduced  matrix  elements  that  is  employed  is  the  matrix 
elements  of  spherical  harmonics.  The  angular  portion  of  our  perturbing  electron-radiator 
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interaction  is  written  in  terms  of  normalized  spherical  harmonics,  defined  as, 
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Several  relations  pertaining  to  reduced  matrix  elements  of  spherical  harmonics  are  writ- 
ten below: 
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6-j  symbols  are  also  used  to  evaluate  the  angular  portions  of  the  matrix  elements. 
The  symbols  are  defined  as  follows: 
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With  regard  to  symmetry  properties,  a 6-j  symbol  is  unchanged  when  any  two  columns 
are  interchanged  and  if  any  two  numbers  in  the  bottom  row  are  interchanged  with  the 
corresponding  two  numbers  in  the  top  row.  The  arguments  in  the  top  row  must  satisfy 
a triangle  relation, 


Ih-jzl  <jl  <72  +73,  |ji  ~ 7.3 1 < J2  < ji  + 7.3,  | Ji  - J2I  < j-i  < Ji  + J2,  (D.15) 
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and  their  sum  must  be  a positive  integer  for  the  6-j  symbol  to  be  nonzero.  Also,  any  set 
of  arguments  that  can  be  permuted  into  the  top  row  must  also  satisfy  these  conditions. 
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